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1 STATE OF THE ART IN LARGE-SCALE SOIL MOISTURE MONITORING 

2 

3 ABSTRACT 

4 Soil moisture is an essential climate variable influencing land — atmosphere interactions, 

5 an essential hydrologic variable impacting rainfall — runoff processes, an essential ecological 

6 variable regulating net ecosystem exchange, and an essential agricultural variable constraining 

7 food security. Large-scale soil moisture monitoring has advanced in recent years creating 

8 opportunities to transform scientific understanding of soil moisture and related processes. These 

9 advances are being driven by researchers from a broad range of disciplines, but this complicates 

10 collaboration and communication. And, for some applications, the science required to utilize 

11 large-scale soil moisture data is poorly developed. In this review, we describe the state of the art 

12 in large-scale soil moisture monitoring and identify some critical needs for research to optimize 

13 the use of increasingly available soil moisture data. We review representative examples of 1) 

14 emerging in situ and proximal sensing techniques, 2) dedicated soil moisture remote sensing 

15 missions, 3) soil moisture monitoring networks, and 4) applications of large-scale soil moisture 

16 measurements. Significant near-term progress seems possible in the use of large-scale soil 

17 moisture data for drought monitoring. Assimilation of soil moisture data for meteorological or 

18 hydrologic forecasting also shows promise, but significant challenges related to model structures 

19 and model errors remain. Little progress has been made yet in the use of large-scale soil 

20 moisture observations within the context of ecological or agricultural modeling. Opportunities 

21 abound to advance the science and practice of large-scale soil moisture monitoring for the sake 

22 of improved Earth system monitoring, modeling, and forecasting. 

23 
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The science and practice of large-scale soil moisture monitoring has entered a stage of 
unprecedented growth with the potential to transform scientific understanding of the patterns and 
dynamics of soil moisture and soil moisture-related processes. Large-scale soil moisture 
monitoring may lead to improved understanding of soil moisture controls on water, energy, and 
carbon fluxes between the land and atmosphere, resulting in improved meteorological forecasts 
and climate projections. Soil moisture measurements are also key in assessing flooding and 
monitoring drought. Knowledge gained from large-scale soil moisture observations can help 
mitigate these natural hazards, yielding potentially great economic and societal benefits. Here 
large-scale refers to spatial support scales of >1 m for an in situ sensor or spatial extents of 
>100 2 km 2 for an in situ network (Crow et al., 2012; Western and Bloschl, 1999). New 
developments continue within the realm of in situ sensors which monitor soil moisture at the 
point-scale, i.e. <1 m support. These point-scale sensors have been reviewed recently 
(Dobriyal et ah, 2012; Robinson et al., 2008) and will not be considered here except within the 
context of large-scale networks. Rather, this review aims to broadly describe the state of the art 
in large-scale soil moisture monitoring. Airborne and satellite remote sensing approaches for 
soil moisture are also considered large-scale monitoring techniques in this review. 

To provide context, it is helpful to begin with a brief historical overview of soil moisture 
monitoring in general. The first major technological advance in modem soil moisture 
monitoring can be traced to the development of the neutron probe after World War II (Evett, 
2001). The measurement of soil moisture based on neutron thermalization first appeared in peer- 
reviewed literature in a paper by Iowa State College (now University) soil physicists, Gardner 
and Kirkham (1952). This technology was soon commercialized under a contract between the US 
Army Corps of Engineers and Nuclear-Chicago Corporation, and by 1960 hundreds of neutron 
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48 probes were in use around the world (Evett, 2001). The neutron probe remained the de facto 

49 standard for indirect soil moisture measurement until a soil physicist and two geophysicists 

50 working for the Government of Canada made a key breakthrough in using dielectric properties to 

51 measure soil water (Topp et al., 1980). Despite initial skepticism from the soil science and 

52 remote sensing communities (Topp, 2006), the time domain reflectometry (TDR) approach of 

53 Topp et al. (1980) eventually became a dominant technology for soil moisture monitoring, and 

54 created for the first time, the possibility of automated, multiplexed, unattended, in situ 

55 monitoring (Baker and Allmaras, 1990). By the 1990s, the TDR technology had proven the 

56 value of electromagnetic methods for monitoring soil moisture, and an avalanche of impedance 

57 or capacitance type probes followed (Robinson et al., 2008). These capacitance probes typically 

58 operate at frequencies much lower than the effective frequency of TDR. As a result these probes 

59 are simpler and less expensive, but also less accurate than TDR (Blonquist et al., 2005). Much 

60 effort has also been devoted to the development of heat dissipation (Fredlund and Wong, 1989; 

61 Phene et al., 1971; Reece, 1996) and heat pulse sensors (Bristow et al., 1993; Campbell et al., 

62 1991; Heitman et al., 2003; Ochsner et al., 2003; Song et al., 1999; Tarara and Ham, 1997) for 

63 soil moisture measurement with reasonable success. 

64 While Canadian researchers were beginning to develop the groundbreaking TDR method, 

65 scientists in the US were pioneering remote sensing of soil moisture from tower, aircraft, and 

66 satellite platforms using microwave radiometers (Schmugge et al., 1974), scatterometers (Dickey 

67 et al., 1974), synthetic aperture radar (Chang et al., 1980), and combined radar/radiometer 

68 systems (Ulaby et al., 1983). A variety of other techniques were also introduced during the same 

69 time, including methods based on polarized visible light (Curran, 1978), thermal inertia (Pratt 

70 and Ellyett, 1979), and terrestrial gamma radiation (Carroll, 1981). Satellite remote sensing 
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approaches in particular have engendered much enthusiasm and interest with their promise of 
global data coverage leading Vinnikov et al. (1999) to speculate that, in regards to long-term soil 
moisture monitoring, “The future obviously belongs to remote sensing of soil moisture from 
satellites”. And, in fact, the intervening decades of research on remote sensing of soil moisture 
are now beginning to bear fruit in terms of operational satellites for large-scale soil moisture 
monitoring. 

Not everyone has been content to wait for the arrival of operational soil moisture 
satellites; rather, some have envisioned and created large-scale in situ monitoring networks for 
soil moisture. The earliest organized networks were in the Soviet Union and used repeated 
gravimetric sampling (Robock et al., 2000). The Illinois Climate Network was the first large- 
scale network to use a nondestructive measurement device, the neutron probe (Hollinger and 
Isard, 1994), while the USDA Natural Resources Conservation Service Soil Climate Analysis 
Network (Schaefer et al., 2007) and the Oklahoma Mesonet (McPherson et al., 2007) pioneered 
the use of automated, unattended sensors in large-scale soil moisture networks during the 1990s. 
Since then numerous networks have emerged around the world, and have come to play vital roles 
in the science and practice of large-scale soil moisture monitoring, not the least of which is their 
role in calibrating and validating satellite remote sensing techniques. 

The past ten years have witnessed the emergence of potentially transformative new soil 
moisture technologies which are beginning to fundamentally alter the possibilities for large-scale 
monitoring. These new methods include the cosmic-ray soil moisture observing system 
(COSMOS), global positioning system (GPS) based techniques, and fiber optic distributed 
temperature sensing (DTS) approaches (Larson et al., 2008; Sayde et al., 2010; Steele-Dunne et 
al., 2010; Zreda et al., 2008). Meanwhile, the number and scope of large-scale automated soil 
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94 moisture monitoring networks has been steadily increasing, both in the US and around the world. 

95 And, in 2009, the European Space Agency launched the Soil Moisture Ocean Salinity (SMOS) 

96 satellite, the first one designed specifically for soil moisture monitoring. 

97 Despite these developments, many challenges remain within the realm of large-scale soil 

98 moisture monitoring. The recent progress in this field has been enabled by contributions from 

99 many different disciplines, and future progress will likely be interdisciplinary, as well. But, 

100 staying informed about new developments can be challenging when the research is spread across 

101 a broad range of science disciplines from soil science to remote sensing to geodesy to 

102 meteorology. Contemporary soil physicists, whose predecessors were instrumental in birthing 

103 the modem era of soil moisture monitoring, have been largely focused on development and 

104 testing of point-scale measurement techniques and have perhaps not been adequately engaged in 

105 advancing the science of large-scale monitoring. Most importantly, the basic science and 

106 technology required to actually use large-scale soil moisture data is relatively under-developed. 

107 There has been a dearth of research investment in developing modeling and forecasting tools 

108 driven by large-scale soil moisture data, especially data from in situ networks. This was 

109 understandable in previous decades when the widespread availability of such data was a distant 

110 prospect, but the circumstances have changed. Soil moisture data are now common and may be 

111 ubiquitous in the near future. 

112 In light of these circumstances, we seek to meet the need for a cross-disciplinary state of 

113 the art review for the sake of improving communication and collaboration. We further seek to 

114 engage and mobilize the expertise of the international soil science, and specifically soil physics, 

115 community in advancing the science and practice of large-scale soil moisture monitoring. Also, 

116 we seek to highlight the pressing need to accelerate the pace of progress in the area of using 
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large-scale soil moisture observations for advanced Earth systems monitoring, modeling, and 
forecasting applications. Our objectives are 1) to succinctly review the state of the art in large- 
scale soil moisture monitoring and 2) to identify some critical needs for research to optimize the 
use of increasingly available soil moisture data. 

This review does not aim to be comprehensive. Rather we have selected specific topics 
which are illustrative of the opportunities and challenges ahead. This review is organized in four 
primary sections: 1) emerging in situ and proximal sensing techniques, 2) dedicated soil moisture 
remote sensing missions, 3) soil moisture monitoring networks, and 4) applications of large-scale 
soil moisture measurements. Some observations regarding primary challenges and opportunities 
for large-scale soil moisture monitoring are provided at the end of the review. 

EMERGING IN SITU AND PROXIMAL SENSING TECHNIQUES 
Soil Moisture Monitoring Using Cosmic Ray Neutrons 

Area-average soil moisture can be measured in the field using cosmic-ray neutron 
background radiation whose intensity in air above the land surface depends primarily on soil 
moisture. The cosmic-ray probe integrates soil moisture over an area hundreds of meters in 
diameter, something that would require an entire network of point measurement devices. 
Measurements can be made using stationary probes, which provide time series of soil moisture, 
or mobile probes, which provide snapshots in time over an area or along a line. 

Cosmic-ray protons that impinge on the top of the atmosphere create secondary neutrons 
that in turn produce additional neutrons, thus forming a self-propagating nucleonic cascade 
(Simpson, 2000; Desilets and Zreda, 2001). As the secondary neutrons travel through the 
atmosphere and then through the top few meters of the biosphere, hydrosphere and lithosphere, 
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fast neutrons are created (Desilets et al., 2010). Because fast neutrons are strongly moderated by 
hydrogen present in the environment (Zreda et al., 2008, 2012), their measured intensities reflect 
variations in the soil moisture (Zreda et al., 2008) and other hydrogen present at and near the 
Earth’s surface (Zreda et al., 2012; Franz et al., 2013). 

The process of neutron moderation depends on three factors that together define the 
neutron stopping power of a material (Zreda et al., 2012): (1) the elemental scattering cross 
section or probability of scattering; hydrogen has a high probability of scattering a neutron; (2) 
the logarithmic decrement of energy per collision, which characterizes how efficient each 
collision is; hydrogen is by far the most efficient element; and (3) the number of atoms of an 
element per unit mass of material, which is proportional to the concentration of the element and 
to the inverse of its mass number. Because of the abundance of water in soils and hydrogen's 
low atomic mass, hydrogen, next to oxygen and silicon, makes up a significant fraction of all 
atoms in many soils. The extraordinarily high stopping power of hydrogen makes the cosmic- 
ray soil moisture method work. 

The fast neutrons that are produced in air and soil travel in all directions within and 
between air and soil and in this way an equilibrium concentration of neutrons is established. The 
equilibrium is shifted in response to changes in the hydrogen content of the media, which in 
practice means changes in the amount of water on or in the soil. Adding water to soil results in 
more efficient moderation of neutrons by the soil, causing a decrease of fast neutron intensity 
above the soil surface. Removing water from the soil has the opposite effect. Thus, by measuring 
the fast neutron intensity in the air the moisture content of the soil can be inferred, for example 
using the equation of Desilets et al. (2010): 


0 = 


a n 


(N/Nj-a, 
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which is plotted in Fig. 1 . In the equation 0 is the neutron-derived moisture content, N is the 
measured neutron intensity, No is the neutron intensity in air above a dry soil (this is a calibration 
parameter obtained from independent in situ soil moisture data), and ao, a i, and a 2 are fitted 
constants that define the shape of the calibration function. Neutron transport modeling shows that 
the shape of the calibration function is similar for different chemical compositions of soil and 
soil textures (Zreda et al., 2008; Desilets et al., 2010) and in presence of hydrogen pools other 
than pore water, for example vegetation or water vapor (Franz et al., 2013; Rosolem et al., 2012). 
Therefore, the same function can be used under different field conditions once corrections are 
made for all pools of hydrogen (Franz et al., 2013). 

The probe senses all hydrogen present within the distance that fast neutrons can travel in 
soils, water, air and other materials near the land surface. That distance varies with the chemical 
composition and density of the material, from centimeters in water through decimeters in soils to 
hectometers in air. The support volume can be visualized as a hemisphere above the soil surface 
placed on top of a cylinder in the soil (Fig. 2). For soil moisture measurements the diameter and 
height of the cylinder are important. The horizontal footprint, which is defined as the area 
around the probe from which 86% (l-e‘ 2 ) of counted neutrons arise, is a circle with a diameter of 
660 m at sea level (Zreda et al., 2008). It decreases slightly with increasing soil moisture content 
and with increasing atmospheric water vapor content, and it increases with decreasing air density 
(decreasing atmospheric pressure or increasing altitude)(Zreda et al., 2012). The horizontal 
footprint has been verified by field measurements (Zweck et al., 201 1). 

The effective depth of measurement, which is defined as the thickness of soil from which 
86% (1-e 2 ) of counted neutrons arise, depends strongly on soil moisture (Zreda et al., 2008). It 
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186 decreases non-linearly from about 70 cm in soils with no water to about 12 cm in saturated soils 

187 and is independent of air density. The effective depth of measurement decreases with increasing 

188 amount of hydrogen in other reservoirs, such as lattice water, soil organic matter or vegetation. 

189 The decrease in the vertical support volume is more significant at the dry end (on the order of 10 

190 cm) than at the wet end (on the order of 1 cm). The vertical footprint has not been verified 

191 empirically. 

192 Neutrons react with any hydrogen present near the Earth’s surface. Therefore, the 

193 measured neutron intensity reflects the total reservoir of neutrons present within the sensing 

194 distance of the probe (Fig. 2), and hence the probe can be viewed as the total surface moisture 

195 probe. The greater the concentration of hydrogen, the greater is its impact on the neutron 

196 intensity. Large near-surface reservoirs of hydrogen, roughly in order of decreasing size, are: (1) 

197 surface water (including snow), (2) soils, (3) lattice water and water in soil organic matter; (4) 

198 vegetation, and (5) atmospheric water vapor. Because the neutron signal integrates all these 

199 factors, isolation of one of these components, for example soil moisture, requires that the others 

200 be: (a) constant in time, (b) if not constant, assessed independently, or (c) negligibly small. In 

201 addition, the support volume (or the measurement volume) will be affected by these other 

202 sources of hydrogen. 

203 Calibration requires simultaneous measurements of area-average soil moisture (0) and 

204 neutron intensity (N), and solving Eq. [1] for the calibration parameter No. Area-average soil 

205 moisture representative of the cosmic-ray footprint is obtained by collecting numerous soil 

206 samples around the cosmic-ray probe and measuring moisture content by the oven-drying 

207 method (Zreda et al., 2012); other methods, such as time-domain reflectometry, can be used as 

208 well. The measured neutron intensities must be corrected for atmospheric water vapor and 
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pressure variations. Soil samples must be analyzed for chemical composition to correct the 
calibration function for any additional water in mineral grains (lattice water) and in organic 
matter present in the soil (Zreda et al., 2012). The presence of that extra water shifts the position 
of the calibration point to the left on the calibration function (Fig. 1), which results in steeper 
curve and thus in reduced sensitivity of neutrons to changes in soil moisture. Other sources of 
water have a similar effect on the calibration function. 

Measurement precision of soil moisture determination is due to neutron counting 
statistics. The counts follow the Poisson distribution (Knoll, 2000) in which for the total number 
of counts, N, the standard deviation is TV 0 ' 5 . Thus, more counts produce better precision (i.e. 
lower coefficient of variation), provided that the neutron intensity remains stationary over the 
counting time. High counting rates are expected under these conditions: (1) high altitude and 
high latitude, because the incoming cosmic-ray intensity, which is the precursor to fast neutrons, 
increases with both (Desilets and Zreda, 2003; Desilets et al., 2006); (2) dry soil, because of the 
inverse relation between soil moisture and neutron intensity (Fig. 1); (3) dry atmosphere, because 
of the inverse relation between atmospheric moisture and neutron intensity (Rosolem et al., 
2012); (4) no vegetation; (5) low lattice and organic matter content of soil. Opposite conditions 
will result in lower counting rates and poorer precision. 

The accuracy of soil moisture determination depends on a few factors related to 
calibration and the presence of other pools of hydrogen within cosmic-ray probe support volume. 
The calibration uncertainty is due to two factors: (1) the accuracy of the independent measure of 

O 'J 

area-average soil moisture, which is usually below 0.01 m m' ; (2) the accuracy of neutron 
count rate at the time of calibration, which is usually around 2%. (These calibration data sets can 
be viewed at cosmos.hwr.arizona.edu.) If these were the only contributing factors, the accuracy 
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232 would be better than 0.01 m 3 m~ 3 . But there are a few complicating factors that may lead to an 

233 increase of the uncertainty. They include atmospheric water vapor, infiltration fronts, changing 

234 horizontal correlation scale of soil moisture, variable vegetation, and variations in the incoming 

235 cosmic-ray intensity. Corrections have been developed for these factors, but their contributions 

236 to the overall uncertainty of soil moisture determination have not been assessed rigorously. 

237 Cosmic-ray soil moisture probes are used as stationary or roving devices. Stationary 

238 probes are installed above the land surface to measure and transmit neutron intensity and 

239 ancillary data at user-prescribed time intervals (Zreda et al., 2012). These measurements are then 

240 used, together with cosmic-ray background intensity data, to compute soil moisture. A network 

241 of stationary probes, called the COsmic-ray Soil Moisture Observing System (COSMOS), is 

242 being installed in the USA, with the main aim to provide area- average soil moisture data for 

243 atmospheric applications (Zreda et al., 2012). Data are available with one hour latency at 

244 http://cosmos.hwr.arizona.edu. Other networks or individual probes are being installed in 

245 Australia (the network named CosmOz), Germany (Rivera Villarreyes et al., 201 1) and 

246 elsewhere around the globe. 

247 A mobile version of the cosmic-ray soil moisture probe, called COSMOS rover, is under 

248 development. Its main application is mapping soil moisture over large areas from a car or an 

249 aircraft; a backpack version is possible as well. The vehicle-mounted instrument is 

250 approximately ten times larger than the stationary cosmic-ray probe to provide more counts 

251 (better statistics) in short time as the vehicle progresses along the route. The measured neutron 

252 intensity is converted to soil moisture using the usual calibration equation (Desilets et al., 2010). 

253 Transects (Desilets et al., 2010) or maps (Zreda et al., 2011) of soil moisture can be produced 

254 within hours or days. Such maps are useful for many applications, including calibration and 
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validation of satellite soil moisture missions SMOS (Kerr et al., 2010) and SMAP (Entekhabi et 
al., 2010). 

Soil Moisture Monitoring Using Global Positioning System Signals 

While the cosmic ray probe utilizes an existing natural “signal”, the ambient fast neutron 
intensity, to infer soil moisture, new methods employing global positioning system (GPS) 
receivers utilize existing anthropogenic signals. The GPS signals follow two types of paths 
between the satellites that transmit GPS signals and the antennas that receive them (Fig. 3). 

Some portion of GPS signals travel directly from satellites to antennas. These direct signals are 
optimal for navigation and geodetic purposes. Antennas also receive GPS signals that reflect off 
the land surface , referred to as multipath by the geodetic community (Georgiadou and 
Kleusberg, 1988). GPS satellites transmit microwave (L-band) signals that are optimal for 
sensing water in the environment (Entekhabi et al., 2010). For bare soil conditions, the reflection 
coefficients depend on permittivity of the soil, surface roughness, and elevation angle of the 
reflections. Therefore, reflected GPS signals can be used to estimate soil moisture, as well as 
other environmental parameters. GPS antennas and receivers can also be mounted on satellites 
(Lowe et al., 2002) or on planes (Katzberg et al., 2005). The data collected by these instruments 
are considered remote sensing observations. Alternatively, GPS reflections can also be measured 
using antennas mounted fairly close to the land surface (Larson et al., 2008; Rodriguez-Alvarez 
et al., 2011a), yielding a hybrid remote sensing-/?? situ observation. Ground-based GPS studies 
use the interference of the direct and reflected GPS signals, and thus the method is often called 
GPS interferometric reflectometry (GPS-IR). 

For GPS-IR systems, the sensing footprint depends on (1) the height of the antenna above 
the ground and (2) the range of satellite elevation angles used in the analysis. As satellite 
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278 elevation angle (E) increases, the portion of the ground that yields specular (i.e. mirror-like) 

279 reflections both shrinks and moves closer to the antenna. For the case of a typical geodetic 

280 antenna height of 2 m, the center of the area sensed varies from 25 m at E= 5 degrees to 5 m at E 

281 = 30 degrees. Larger sampling areas can be achieved by raising the antenna to heights of -100 

282 m, above which observations are complicated by the GPS code lengths (Rodriguez- Alvarez et 

283 al., 2011a). As GPS is a constellation of more than 30 satellites, different GPS satellites rise and 

284 set above a GPS soil moisture site throughout the day. These reflections are measured from 

285 different azimuths depending on the orbital characteristics of each satellite. For the best sites, 

286 more than 60 soil moisture estimates can be made per day. So, the soil moisture data estimated 

287 from GPS reflections should be considered as daily in temporal frequency, once averaged over 

288 an area of -1000 m for antenna heights of 2 m (Larson et ah, 2008). 

289 Two methods of GPS soil moisture sensing are currently being developed. The first is 

290 based on using GPS instruments designed for geodesists and surveyors. These GPS instruments 

291 traditionally measure the distance between the satellites and antenna in order to estimate 

292 position. However these GPS instruments also measure signal power, or signal-to-noise ratio 

293 (SNR). Embedded on the direct signal effect are interference fringes caused by the reflected 

294 signal being in or out of phase with respect to the direct signal. The SNR frequency is primarily 

295 driven by the height of the antenna above the ground. As permittivity of the soil changes, the 

296 amplitude, phase, and frequency of the SNR interferogram varies. (Larson et al., 2010; 

297 Zavorotny et al., 2010). Of the three parameters, the phase of the SNR interferogram is the most 

298 useful for estimating soil moisture. 

299 Chew et al (2013) demonstrated theoretically that phase varies linearly with surface soil 

300 moisture. For the soils described by Hallikainen et al. (2005), the slope of this relationship does 
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not vary with soil type. For most conditions, phase provides a good estimate of average soil 
moisture in the top 5 cm. The exception is when very wet soil overlies dry soil, for example 
immediately following short-duration rainstorms when the wetting front has not propagated to ~5 
cm (Larson et al., 2010). Estimates of soil moisture from phase have been compared to in situ 
soil moisture measurements (Fig. 4). At grass-dominated sites with relatively low vegetation 
water content (<0.5 kg m' ), SNR phase varies linearly with in situ soil moisture (r >0.8) (Larson 
et al., 2008; and unpublished data) as predicted by Chew et al. (2013). The vegetation at these 
sites is typical of many rangeland areas in the western U.S. A SNR interferogram is also 
affected by higher water content vegetation, for example that which exists in irrigated 
agricultural fields (Small et al., 2010). Methods are being developed to retrieve surface soil 
moisture from SNR interferograms under these conditions. 

One advantage to using geodetic GPS equipment to measure soil moisture is that existing 
geodetic networks can provide much needed hydrologic information. The National Science 
Foundation’s Plate Boundary Observatory (PBO) network has more than 1 100 stations with 
effectively identical GPS instrumentation. Many of the stations are located amidst complex 
topography, which does not facilitate estimation of soil moisture via GPS-IR. However, soil 
moisture is being estimated at 59 stations with relatively simple topography. The data is updated 
daily and is available at http://xenon.colorado.edu/portal/. 

A second GPS soil moisture sensing method is also under development (Rodriguez- 
Alvarez et al. 2009). Similar to Larson et al. (2008), this system measures the interference 
pattern resulting from the combination of direct and reflected GPS signals. A dual polarization 
antenna measures power of the vertically- and horizontally-polarized signals separately, which is 
not possible using standard geodetic instrumentation. The satellite elevation angle at which 
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324 reflectivity of the vertically-polarized signal approaches zero, i.e. the Brewster angle, varies with 

325 soil moisture (Rodriguez- Alvarez et al., 2011a). The existence of this Brewster angle yields a 

326 notch in the interference pattern. The position of the notch is then used to infer soil moisture. 

327 Over a bare soil field, this technique yielded 10 soil moisture estimates over a one month 

328 period; they show good agreement with those measured in situ at a depth of 5 cm (RMS error < 

329 0.03) (Rodriguez-Alvarez et al.,2009). A vegetation canopy introduces additional notches to the 

330 observed interference pattern. The position and amplitude of these notches can be used to infer 

331 both vegetation height and soil moisture. This approach yielded excellent estimates of com 

332 height throughout a growing season (RMS error = 6.3 cm) (Rodriguez-Alvarez et al., 2011b). 

333 Even beneath a 3m tall com canopy, soil moisture estimates typically differed by <0.04 m m' 

334 from those measured with in situ probes at 5 cm. The main difference between these two 

335 approaches is that the approach of Larson et al. (2008) uses commercially-available geodetic 

336 instrumentation - which typically already exists - and can be simultaneously used to measure 

337 position. The approach of Rodriguez-Alvarez et al. (2009) uses a system specifically designed 

338 for environmental sensing, but it is not yet commercially-available, 

339 Soil Moisture Monitoring Using Distributed Temperature Sensing 

340 Much as the Larson et al. (2008) GPS-IR method repurposes commercially available GPS 

341 receivers to monitor soil moisture, other researchers have sought to develop new soil moisture 

342 monitoring methods using commercially available distributed temperature sensing (DTS) 

343 systems. In a DTS system, an optical instrument is used to observe temperature along a 

344 continuum of points within an attached optical fiber cable, typically by the principle of Raman 

345 scattering (Selker et al., 2006). The spatial location corresponding to each temperature 

346 measurement is determined based on the travel time of light in the fiber in a manner analogous to 
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TDR. Weiss (2003) pioneered the use of DTS systems for soil moisture monitoring by 
successfully demonstrating the potential use of fiber optics to detect the presence of moisture in a 
landfill cover constructed from sandy loam soil. A 120-V generator supplied current to the 
stainless steel sheath of a buried optical fiber cable for -626 s at a rate of 18.7 W m' 1 , and the 
corresponding spatially variable temperature rise of the cable was observed at 40-s temporal 
resolution and 1-m spatial resolution. Analysis of the temperature rise data using the single 
probe method (Carslaw and Jaeger, 1959) resulted in satisfactory estimates of the spatial 
variability of soil thermal conductivity along the cable, which in turn reflected the imposed 
spatial variability of soil moisture. However, the temperature uncertainty achieved was ~0.55°C, 
and Weiss concluded that without improvements in signal- to-noise ratio, that system would not 
be able to resolve small changes in soil moisture above 0.06 m 3 m~ 3 for the sandy loam soil used 
in that study. 

The potential of using passive (unheated) DTS methods for soil moisture estimation was 
explored by Steele-Dunne et al. (2010). Optical fiber cable was installed in a tube on the soil 
surface and at depths of 8 and 10 cm. The soil texture was loamy sand, and the vegetation cover 
was sparse grass. With temperatures from the upper and lower cables as time-dependent 
boundary conditions, the temperature at the middle cable was modeled by numerical solution of 
the 1-D heat conduction equation. A numerical search routine was used to find the thermal 
diffusivity which produced the best agreement between the simulated and observed temperatures 
at the 8 cm depth. The results demonstrated that the passive DTS system could detect temporal 
changes in thermal diffusivity associated with rainfall events, but the accuracy of the diffusivity 
estimates was hindered by uncertainties about the exact cable depths and spacings. Furthermore, 
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369 deriving soil moisture estimates was complicated by uncertainty and nonuniqueness in the 

370 diffusivity — soil moisture relationship. 

371 Sayde et al. (2010) modified the active DTS approach of Weiss (2003) by interpreting the 

372 temperature rise data in terms of cumulative temperature increase, i.e. the integral of the 

373 temperature rise from the beginning of heating to some specified time limit. Based on a 

374 laboratory sand column experiment with 2-min, 20 W m' 1 heat pulses, they developed an 

375 empirical calibration function which fit the observed cumulative temperature increase (0 to 120 

376 s) versus soil moisture data. Based on that function and the observed uncertainty in the 

377 cumulative temperature increase data, the uncertainty in the soil moisture estimates would 

O O ... O 

378 increase approximately linearly from 0.001 m m' when soil moisture is 0.05 m m' to 0.046 m 

'J 'J 'J 

379 m' when soil moisture is 0.41 m m' . Gil-Rodriguez et al. (2012) used the approach of Sayde 

380 et al. (2010) to satisfactorily monitor the dimensions and evolution of the wetted bulb during 

381 infiltration beneath a drip emitter in a laboratory column of sandy loam soil. 

382 Striegl and Loheide (2012) used an active DTS approach to monitor spatial and temporal 

383 dynamics of soil moisture along a 130-m transect associated with a wetland reconstruction 

384 project (Fig. 5). They used a 10-min, 3 W m' 1 heat pulse, a lower heating rate than used in 

385 previous active DTS studies. They followed Sayde et al. (2010) in adopting a primarily 

386 empirical calibration approach, but rather than cumulative temperature increase, they related soil 

387 moisture to the average temperature rise observed from 380 to 580 s after the onset of heating. A 

388 calibration function was developed by relating the observed temperature rise data to independent 

389 soil moisture measurements at three points along the transect, and the resulting function had a 

390 RMSE = 0.016 m 3 m' 3 for soil moisture < 0.3 1 m 3 m' 3 but RMSE = 0.05 m 3 m' 3 for wetter 

391 conditions. Their system successfully monitored field scale spatiotemporal dynamics of soil 
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moisture at 2-m and 4-h resolution across a 2-mo period consisting of marked wetting and drying 
cycles (Fig. 6). 

The passive and active DTS methods for monitoring soil moisture offer the potential for 
unmatched spatial resolution (<1 m) in long-term soil moisture monitoring on field scale (>100 
m) transects. These methods may in the near future greatly impact our understanding of the fine- 
scale spatiotemporal structure of soil moisture and shed new light on the factors influencing that 
structure. Thus far, the active DTS methods have shown more promise than passive DTS, but 
more sophisticated data assimilation approaches for interpreting passive DTS data are in 
development. The active DTS method is still in its infancy, and many key issues remain to be 
addressed. None of the active DTS methods developed to date involve spatial variability in the 
soil moisture calibration function, so heterogeneity in soil texture and bulk density could give 
rise to appreciable uncertainties in field settings. Field installation of the optical fiber cables at 
the desired depths with good soil contact and minimal soil disturbance is also a significant 
challenge. Custom designed cable plows (Steele-Dunne et al., 2010) and commercial vibratory 
plows (Striegl and Loheide, 2012) have been used with some success. The active DTS methods 
have demonstrated good precision for low to moderate soil moisture levels, but further 
improvements in measurement precision are needed for wet conditions. Obtaining good quality 
temperature measurements using a DTS instrument in the field requires that thermally-stable 
calibration baths be included in the system design. The instrument itself must also be in a 
thermally-stable environment because sizeable errors can result from sudden changes in the 
instrument temperature (Striegl and Loheide, 2012). The measurement principles behind DTS 
are discussed in more detail by Selker et al. (2006), and practical aspects of DTS, including key 
limitations and uncertainties, are described by Tyler et al. (2009). 
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415 

416 DEDICATED SOIL MOISTURE REMOTE SENSING MISSIONS 

417 Soil Moisture and Ocean Salinity (SMOS) 

418 The Soil Moisture and Ocean Salinity mission (Kerr et al., 2010), an Earth Explorer 

419 Opportunity mission, was successfully launched on November 2, 2009, and successfully 

420 concluded its commissioning phase in May 2010. It was developed under the leadership of the 

421 European Space Agency (ESA) with the Centre National d’ Etudes Spatiales (CNES) in France 

422 and the Centro para el Desarrollo Teccnologico Industrial (CDTI) in Spain. 

423 Microwave radiometry at low frequencies is an established technique for estimating 

424 surface soil moisture with an adequate sensitivity. The choice of L-band as the spectral range in 

425 which to operate was determined from a large number of studies that demonstrated L-band has 

426 high sensitivity to changes of moisture in the soil (Schmugge and Jackson, 1994) and salinity in 

427 the ocean (Lagerloef, 2001). Furthermore, observations at L-band are less susceptible to 

428 attenuation due to the atmosphere or the vegetation than measurements at higher frequencies 

429 (Jackson and Schmugge, 1989; Jackson and Schmugge, 1991). L-band also enables a larger 

430 penetration depth into the surface soil layer than is possible with shorter wavelengths 

431 (Escorihuela et al., 2010) 

432 Even though the L-band radiometry concept was demonstrated early by a space 

433 experiment (SKYLAB) back in the 1970’s, no dedicated space mission followed because 

434 achieving a suitable ground resolution (< 50-60 km) required a prohibitive antenna size (> 4 m). 

435 The so-called interferometry design, inspired from the very large baseline antenna concept (radio 

436 astronomy), made such a venture possible. The idea consists of deploying small receivers in 

437 space (located on a deployable structure), then reconstructing a brightness temperature (T B ) field 
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with a resolution corresponding to the spacing between the outmost receivers. The two- 
dimensional interferometer allows measuring T B at several incidence angles, with full 
polarization. Such an instrument instantaneously records a whole scene; as the satellite moves, a 
given point within the 2D field of view is observed from different view angles. The series of 
independent measurements allows retrieving surface parameters with much improved accuracy. 

The baseline SMOS payload is thus an L-band (1.413 GHz, 21 cm - located within the 
protected 1400-1427 MHz band) 2D interferometric radiometer that is Y shaped with three 4.5 m 
arms as shown in Figure 7. SMOS is on a sun synchronous (6 a.m. ascending) circular orbit and 
measures the brightness temperature emitted from the Earth at L-band over a range of incidence 
angles (0 to 55°) across a swath of approximately 1000 km (covering the globe twice in less than 
3 days) with a spatial resolution of 35 to 50 km (average is 43 km)(Kerr et al., 2001; Kerr et al., 
2010 ). 

The SMOS mission originated from recognition of the need for accurate, global, soil 
moisture monitoring from space. Short wave radiation instruments were quickly discarded 
because of poor sensitivity and the negative impact of cloud cover (Kerr, 2007). Use of thermal 
infra-red also suffered complications due to the need for accurate knowledge of forcings (Kerr, 
2007). Radars and synthetic aperture radar (SAR) typically suffer from low temporal resolution, 
often compensated by a high spatial resolution. Another limitation of these active techniques is 
linked to the difficulty in separating the surface roughness contribution from that of soil 
moisture, often requiring the “change detection approach” (Moran et al., 1998; Moran et al., 
2002). Another possibility is to use scatterometers which are characterized by a lower spatial 
resolution but higher temporal resolution adequate for water budget monitoring. The European 
Remote Sensing Satellite 1 (ERS-1), European Remote Sensing Satellite 2 (ERS-2), and then 
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461 MetOp scatterometers offered such opportunities (Magagi and Kerr, 1997a; Magagi and Kerr, 

462 1997b; Magagi and Kerr, 2001; Wagner et al., 2007) relying on a change detection approach, and 

463 thus not delivering absolute values. Consequently, it seemed logical to investigate passive 

464 microwaves at low frequencies as the ultimate approach to infer soil moisture from space with 

465 the caveat of lower spatial resolution. Interferometry was first put forward by D. LeVine et al. in 

466 the 1980’s (the ESTAR project) and validated with an airborne prototype (Le Vine et al., 1994; 

467 Le Vine et al., 1990). In Europe, an improved concept was next proposed to the European Space 

468 Agency (ESA): the Microwave Imaging Radiometer using Aperture Synthesis (MIRAS) concept 

469 (Goutoule, 1995). This concept has now materialized into the SMOS mission. 

470 The SMOS data have demonstrated good sensitivity and stability. The data quality was 

471 sufficient to allow the production - from an interferometer - of prototype global soil moisture 

472 maps within one year after launch. It was the first time ever such maps were obtained. Initially, 

473 the accuracy was relatively poor and many retrievals were not satisfactory. The data were much 

474 impaired by radio frequency interference (RFI) leading to degraded measurements in several 

475 areas including parts of Europe and China (Oliva et al., 2012). With the help of the SMOS team, 

476 ESA and CNES took actions to reduce RFI. Specific RFI sources are now identified and 

477 localized then provided to ESA personnel who interact directly with the appropriate national 

478 agencies. These efforts have resulted in over 215 powerful and persistent RFI sources 

479 disappearing, including the US Defense Early Warning System in Northern Canada and many 

480 sources in Europe. Unfortunately, the remaining number of sources in some countries is large. 

481 While RFI reduction and retrieval algorithm improvement efforts were ongoing, first 

482 attempts to use SMOS data in a variety of applications were investigated. The first topic was to 

483 validate the soil moisture retrievals against in situ measurements, model outputs or other remote 
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sensing platforms. These efforts showed that SMOS soil moisture retrievals equaled or 
surpassed the best techniques previously available with ample room for improvements (A1 Bitar 
et al., 2012; Albergel et al., 2012a; Bircher et al., 2012; Jackson et al., 2012; Kerr et al., 2012; 
Leroux et al., 2012; Mecklenburg et al., 2012; Rahmoune et al., 2012; Schwank et al., 2012). 
Floods in Pakistan occurring just after the end of the commissioning phase proved that SMOS 
was able to track such events in spite of the complex topography. The floods in the US during 
spring 2011 were clearly seen in the SMOS data, as well as the related human activities such as 
levee bursting. Most of the large flood events occurring since launch have been monitored, and 
SMOS has shown its ability to provide information quickly and regularly, not being hindered by 
either cloud cover or revisit time, at the cost of a spatial resolution which is lower than optimal 
for this application. In several cases, the arrival of intensive rains (Yasi Huricane in Australia for 
instance) SMOS data enabled anticipation of flooding risks as a function of soil wetness prior to 
the rains. 

Currently intensive work is underway to improve the spatial resolution of the SMOS 
retrievals with good success using disaggregation techniques (Merlin et al., 2010; Merlin et al., 
2012). Current activities are also devoted to the estimation of water in the entire root zone with 
some success, to inferring a drought index, as well as to the possibilities of using SMOS for 
routing modeling (Pauwels et al., 2012) and for correcting space estimates of rainfall over land. 
Work to address science challenges affecting the SMOS data is also ongoing (Kerr et al., 2012). 
One may cite for instance improving knowledge of water bodies and their temporal evolution, 
modeling of forests, improving knowledge of soil texture on a global basis and - of course - 
general instrument calibration issues. Other current efforts are devoted to improving the auxiliary 


22 

5585 Guilford Rd„ Madison Wl 53711 


Page 23 of 110 


Soil Science Society of America Journal 


506 data sets used in retrievals (e.g. snow and frozen soils) as well as improving underlying models 

507 (e.g. dielectric permittivity, forest emissions, etc. . .). 

508 Currently, SMOS data is freely available from different sources, depending on the type 

509 (or Level) of data required. Level 1 (brightness temperatures) and Level 2 (ocean salinity over 

510 oceans or soil moisture/ vegetation opacity over land) data are available from ESA. The data is 

511 provided in swath mode (half orbits from pole to pole) in BinHex format and on the ISEA 49H 

512 grid. These Levels are available through the ESA (https://earth.esa.int/web/guest/missions/esa- 

513 operational-eo-missions/SMOS). Level 3 data consist of composited data over either one day 

514 (i.e. all the Level 2 data of one day in the same file), three days, ten days, or one month and over 

515 the globe (either morning or afternoon passes) for soil moisture and vegetation opacity. Over 

516 oceans the sampling is either daily or monthly. Level 3 data are available from the Centre Aval 

517 De Traitement des donnees SMOS (CATDS) through an ftp site (ftp://eftp.ifremer.fr/catds/cpdc; 

518 write to support@catds.fr to get access). The data is provided in NetCDF format on the EASE 

519 grid (25km sampling). Other soil moisture products (root zone soil moisture, drought indices, 

520 etc. . .) will soon be available from the same site. Finally, the implementation of these Level 3 

521 products is expected to bring significant improvements, particularly in the vegetation opacity 

522 retrieval using temporal information (Jacquette et al., 2010). Figure 8 shows a typical monthly 

523 Level 3 soil moisture product. 

524 After the successful launch of SMOS, Aquarius was successfully launched on June 10 

525 2011 and SMAP (see below) is scheduled to launch in 2014. These NASA missions are in a way 

526 complementary to SMOS and should also bring their yield of good results. New breakthroughs 

527 are expected either using single instrument measurements or, more likely, through synergisms 

528 with other sensors either in the optical/thermal infra-red range or with active/ passive microwave 
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sensors. But, a lingering challenge remains. How to achieve better temporal and spatial sampling 
of the globe for soil moisture? The simplest approach relies on dis-aggregation techniques. 

These techniques use data from high resolution sensors to distribute soil moisture as measured by 
an interferometer, and successful results have been already obtained (Merlin et al., 2010; Merlin 
et al., 2005; Merlin et al., 2012). Recognizing the challenge of improving spatial resolution, 
CNES has initiated research activities whose goal is to develop a new mission which would 
fulfill all the SMOS requirements but with a ten times better spatial resolution and an improved 
sensitivity (factor of three for salinity applications), paving the way to more applications in water 
resources management, coastal area monitoring, basin hydrology or even thin sea ice monitoring 
(Kaleschke et al., 2012). The concept, named SMOSNEXT, is based of merging spatial and 
temporal 2D interferometry and is currently undergoing phase 0 at CNES with a proof of concept 
experiment funded by the R&D program. 

Soil Moisture Active/Passive Mission (SMAP) 

The NASA Soil Moisture Active Passive (SMAP) mission (Entekhabi et al., 2010) is 
scheduled to launch in October 2014. Like SMOS, the SMAP mission will utilize L-band 
measurements to determine surface soil moisture conditions, but SMAP will feature both active 
and passive L-band instruments, unlike SMOS which relies on passive measurements alone. The 
SMAP measurement objective is to provide frequent, high-resolution global maps of near- 
surface soil moisture and freeze/thaw state. These measurements will greatly improve estimates 
of water, energy and carbon fluxes between the land and atmosphere. Observations of the timing 
of freeze/thaw transitions over boreal latitudes will help reduce major uncertainties in 
quantifying the global carbon balance. The SMAP soil moisture mission requirement is to 
provide estimates of soil moisture at 10 km spatial resolution in the top 5 cm of soil with an error 
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552 of no greater than 0.04 cm 3 cm' 3 at three-day average intervals over the global land area, 

553 excluding regions of snow and ice, frozen ground, mountainous topography, open water, urban 

554 areas, and vegetation with water content greater than 5 kg m' (averaged over the spatial 

555 resolution scale). This level of performance will enable SMAP to meet the needs of 

556 hydrometeorology and hydroclimate applications. 

557 The SMAP spacecraft (Fig. 9) will carry two L-band microwave instruments: a non- 

558 imaging synthetic aperture radar operating at 1.26 GHz and a digital radiometer operating at 1.41 

559 GHz. The instruments share a rotating 6-meter offset-fed mesh reflector antenna that sweeps out 

560 a 1000 km- wide swath. The spacecraft will operate in a 685-km polar orbit with an 8-day 

561 repeating ground track. The instrument is designed to provide high-resolution and high-accuracy 

562 global maps of soil moisture at 10 km resolution and freeze/thaw state at 3 km resolution, every 

563 two to three days using combined active (radar) and passive (radiometer) instruments. The 

564 radiometer incorporates radio-frequency interference (RFI) mitigation features to protect against 

565 RFI from man-made transmitters. The radiometer provides high soil moisture accuracy at 

566 moderate spatial resolutions (40 km) by measuring microwave emission from the surface. The 

567 emission is relatively insensitive to surface roughness and vegetation as compared to the radar. 

568 The radar measures backscatter from the surface with high spatial resolution (1-3 km in high 

569 resolution mode), but is more influenced by roughness and vegetation than the radiometer. The 

570 combined radar and radiometer measurements provide soil moisture accuracy approaching 

571 radiometer-based retrievals but with intermediate spatial resolution approaching radar-based 

572 resolutions. Thus, the driving aspects of SMAP’s measurement requirements include 

573 simultaneous measurement of L-band brightness temperature and backscatter with a three-day 


25 

5585 Guilford Rd„ Madison Wl 53711 


574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 


Soil Science Society of America Journal 


Page 26 of 110 


revisit and high spatial resolution (40 km and 3 km, respectively). The combined SMAP soil 
moisture product will be produced at 10-km resolution. 

The planned data products for SMAP include: Level IB and 1C instrument data 
(calibrated and geolocated radar backscatter cross sections and radiometer brightness 
temperatures); Level 2 geophysical retrievals of soil moisture; Level 3 daily composites of Level 
2 surface soil moisture and freeze/thaw state data; and Level 4 value-added data products that are 
based on assimilation of SMAP data into land surface models. The SMAP Level 1 radar data 
products will be archived and made available to the public by the Alaska Satellite Facility in 
Fairbanks, AK, while the Level 1 radiometer and all higher level products will be made available 
by the National Snow and Ice Data Center in Boulder, CO. 

The Level 4 products will support key SMAP applications and address more directly the 
driving science questions of the SMAP mission. SMAP L-band microwave measurements will 
provide direct sensing of surface soil moisture in the top 5 cm of the soil column. However, 
several of the key applications targeted by SMAP require knowledge of root zone soil moisture 
in the top 1 m of the soil column, which is not directly measured by SMAP. The SMAP Level 4 
data products are designed to fill this gap and provide model-based estimates of root zone soil 
moisture that are informed by and consistent with assimilated SMAP surface observations. Error 
estimates for the Level 4 soil moisture product will be generated as a by-product of the data 
assimilation system. A Level 4 carbon product will also be produced that utilizes daily soil 
moisture and temperature inputs with ancillary land cover classification and vegetation gross 
primary productivity (GPP) inputs to compute the net ecosystem exchange (NEE) of carbon 
dioxide with the atmosphere over northern (> 45°N latitude) vegetated land areas. The NEE of 
carbon dioxide with the atmosphere is a fundamental measure of the balance between carbon 
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597 uptake by vegetation GPP and carbon losses through autotrophic and heterotrophic respiration. 

598 The SMAP Level 4 carbon product will provide regional mapped measures of NEE and 

599 component carbon fluxes that are within the accuracy range of tower-based eddy covariance 

600 measurement approaches. 

601 Airborne Microwave Observatory of Subcanopy and Subsurface Mission (AirMOSS) 

602 Current estimates of NEE at regional and continental scales contain such important 

603 uncertainties that amongst the 11 or so models tested there could be differences of 100 percent or 

604 more, and it is not always clear whether the North American ecosystem is a net sink or source for 

605 carbon (Denning et al., 2005; Friedlingstein et al., 2006). Root zone soil moisture (RZSM) is 

606 widely accepted to have a first-order effect on NEE (Law et al., 2002), yet RZSM measurements 

607 are not often available with spatial or temporal extent necessary for input into regional or 

608 continental scale NEE models. Unlike the L-band missions, SMOS and SMAP, which measure 

609 surface soil moisture, the AirMOSS mission is designed to measure RZSM directly. The 

610 hypothesis of the NASA-funded AirMOSS project is that integrating spatially and temporally 

611 resolved observations of root zone soil moisture into ecosystem dynamics models can 

612 significantly reduce the uncertainty of NEE estimates and carbon balance estimates. 

613 The AirMOSS plan is to provide measurements to estimate RZSM using an ultra-high 

614 frequency (UHF - also referred to as P-band) airborne radar, over representative sites of the nine 

615 major North American biomes (Fig. 10). These include: boreal forest (Biome 1); temperate 

616 grassland and savanna shrublands (Biome 5); temperate broadleaf and mixed forest (Biome 2); 

617 temperate conifer forest east (Biome 3); temperate conifer forest west (Biome 4); Mediterranean 

618 woodlands and shrublands (Biome 6); arid and xeric shrublands (Biome 7); tropical and 

619 subtropical dry forest (Biome 8); and tropical and subtropical moist forest (Biome 9). These 
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radar observations will be used to retrieve root zone soil moisture, which along with other 
ancillary data, such as topography, land cover, and various in-situ flux and soil moisture 
observations, will provide the first comprehensive data set for understanding the processes that 
control regional carbon and water fluxes. The public access web site for the AirMOSS project is 
http://airmoss.jpl.nasa.gov/. 

The airborne P-band radar system, flown on a NASA Gulfstream III aircraft, has a flight 
configuration over the experimental sites of typically 100 km by 25 km made up of four flight 
lines (Fig. 11). This represents an intermediate footprint between the flux tower observations (on 
the order of 1 km) and regional to continental scale model simulations. It should be noted that 
each AirMOSS flux site also has a hydrologic modeling domain of on the order of 100 km by 
100 km that will be populated with the corresponding ancillary data sets to allow flexibility in 
the flight line design. The hydrologic simulation domain is determined based on maximizing the 
overlap of full watersheds with the actual flight domain. These watersheds are to be simulated 
using the fully distributed, physically-based finite element model PIHM (Penn State Integrated 
Hydrologic Model) (Qu and Duffy, 2007; Kumar et al., 2010). Carbon dioxide modeling will be 
performed using the Ecosystem Demography (ED2) model (Moorcroft et al., 2001). Each 
AirMOSS site has flux tower measurements for water vapor and carbon dioxide made using an 
eddy covariance system. 

The P-band radar operates in the 420 to 440 MHz frequency range (70 cm), with a longer 
wavelength than typically used in the L-band missions such as SMOS or the upcoming U.S. 
SMAP mission (next section). Previous studies using similar wavelengths have shown that 
RZSM can be computed with an absolute accuracy of better than 0.05 m m' and relative 

3 3 1 

accuracy of 0.01 to 0.02 m m' through a canopy of up to 120 Mg ha' and to soil depths of 50 to 
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643 100 cm, depending on the vegetation and soil water content (Moghaddam et al., 2000; 

644 Moghaddam 2009). This P-band radar system has evolved from the existing Uninhabited Aerial 

645 Vehicle Synthetic Aperture Radar (UAVSAR) subsystems, including the radio frequency 

646 electronics subsystem (RFES), the digital electronics subsystem (DES), the power subsystem, 

647 and the differential GPS subsystem. In fact, the P-band radar system is mounted within the 

648 UAVSAR platform pod on the NASA Gulfstream III thereby negating the requirement for 

649 additional air-worthiness trials. The radar backscatter coefficients are available at both 0.5 arc- 

650 second (approximately 15 m, close to the fundamental spatial resolution of the radar) and at 3 

651 arc-second (approximately 100 m), and the retrieved root zone soil moisture maps will be at 3 

652 arc-second resolution. 

653 AirMOSS flight operations began in Fall of 2012, and all sites in North America except 

654 the tropical sites (Chamela, Mexico and La Selva, Costa Rica) and the woody Savanna site 

655 (Tonzi Ranch, CA) were flown. These P-band data are currently undergoing initial calibration. 

656 However, a three-band raw data image showing the spatial variation of soil moisture over the 

657 Metolius, Oregon site, along with soil roughness and vegetation effects which have not yet been 

658 removed, is shown in Fig. 12. 

659 

660 LARGE-SCALE SOIL MOISTURE MONITORING NETWORKS 

661 Soil moisture networks with spatial extents of >100 km are well-suited for monitoring 

662 the meteorological scale of soil moisture spatial variability as defined by Vinnikov et al. (1999) 

663 because atmospheric forcings often exhibit spatial autocorrelation lengths of 100s of km. These 

664 large-scale networks are also appropriate for studies related to basin-scale hydrology and meso- 

665 scale meteorology. Numerous smaller networks exist worldwide with spatial extents <100 km , 
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both within and outside the US. For example, the USD A Agricultural Research Service (ARS) 
has developed several soil moisture networks to enhance their experimental watershed program. 
Locations include the Little Washita in Oklahoma, Walnut Gulch in Arizona, Reynolds Creek in 
Idaho, and Little River in Georgia (Jackson et al., 2010). The smaller scale networks are often 
well-suited for watershed-scale hydrologic studies. A recent surge in the creation of these 
smaller-scale networks has been driven by the need to validate soil moisture estimates from 
satellites such as SMOS and SMAP. A partial list of current and planned soil moisture networks 
with spatial extents <100 km was provided by Crow et al. (2012). 

Large-Scale Soil Moisture Networks in the United States 

Large-scale soil moisture networks in the U.S. are currently operating in a variety of 
configurations at both national and state levels (Fig. 13, Table 1). In 1981, the Illinois Water 
Survey began a long term program to monitor soil moisture in situ (Hollinger and Isard, 1994; 
Scott et al., 2010). This network was limited by its use of neutron probes, which required 
significant resources to operate and maintain. These neutron probes were used to measure soil 
moisture as frequently as twice a month. These stations were collocated with the Illinois Climate 
Network stations as the Water and Atmospheric Resources Monitoring Program and ultimately 
totaled 19 stations with measurements from the surface to a depth of 2 m. Beginning in 1998, 
these stations were converted to continuously monitor soil moisture using dielectric sensors 
(Hydra Probe, Stevens Water Monitoring Systems, Inc., Portland, OR), providing regular 
statewide estimates of soil moisture. 

The next network to develop was in Oklahoma, which has become an epicenter of 
mesoscale weather and climate research. The Oklahoma Mesonet was launched in 1991 and 
became fully operational in 1994, now consisting of 120 stations, with at least one station in each 
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689 county of Oklahoma (Brock et ah, 1995; McPherson et ah, 2007). Each station hosts a suite of 

690 meteorological measurements, including air temperature, wind speed and direction, air pressure, 

691 precipitation, and soil temperature. These stations monitor soil matric potential using heat 

692 dissipation sensors (CS-229, Campbell Scientific, Inc., Logan, UT) at the 5 cm, 25 cm, and 60 

693 cm depths, with archived data from the 75 cm depth available for some sites. These matric 

694 potentials can be converted to soil moisture estimates via site- and depth-specific water retention 

695 curves (Illston et ah, 2008). Also distributed through Oklahoma is a network of stations 

696 belonging to the Southern Great Plains (SGP) site of the US DOE Atmospheric Radiation 

697 Measurement Program (Schneider et ah, 2003). This network uses the same type of sensor as the 

698 Oklahoma Mesonet. This network began in 1996 and spans portions of Oklahoma and Kansas. 

699 There are a variety of facilities administered by the ARM-SGP site including a large central 

700 facility, as well as extended and boundary facilities, hosting a variety of meteorological, surface, 

701 and soil profile measurements. 

702 While the Oklahoma Mesonet was being developed, the USDA Natural Resource 

703 Conservation Service (NRCS) began a pilot soil moisture/soil temperature project to monitor 

704 these parameters on a national scale. This project developed into the Soil Climate Analysis 

705 Network (SCAN), which now numbers approximately 180 stations across the U.S. (Schaefer et 

706 ah, 2007). This network has a standardized depth profile of Hydra Probe sensors at 5, 10, 20, 50, 

707 and 100 cm. A similar network to SCAN is the Climate Reference Network (CRN), operated by 

708 the NOAA National Climatic Data Center (Palecki and Groisman, 2011). This network 

709 commissioned 114 stations to provide a national scale weather and climate monitoring network. 

710 Soil moisture sensors are being added to these stations currently based on the SCAN 

711 configuration (Hydra Probes at 5, 10, 20, 50, and 100 cm), but three profiles of sensors are 
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installed at each site providing data in triplicate for each depth. In addition to soil moisture, 
standard weather variables such as air temperature, solar radiation, precipitation, and wind speed 
are also collected. 

A number of other state-wide or large-scale networks have been developed since the mid 
1990s. In 1998, the High Plains Regional Climate Center added soil moisture sensors to 14 
Automated Weather Data Network (AWDN) stations in Nebraska. Since then sensors have been 
added to other stations, so now there are 53 stations throughout the state monitoring soil moisture 
on an hourly basis. These stations monitor soil moisture using impedance sensors (Theta Probe 
ML2x, Delta-T Devices, Ltd., Cambridge, UK) at depths of 10 cm, 25 cm, 50 cm, and 100 cm 
(Hubbard et al., 2009). 

The North Carolina Environment and Climate Observing Network (ECONet) has been in 
operation since 1999 when 27 stations were instrumented with Decagon ECHO probes (Weinan 
et al., 2012). In 2003, these stations were converted to Theta Probe sensors and the network was 
expanded to 37. Unlike most other networks, this network does not have a near-surface 
measurement depth as these data are collected only at a 20 cm depth. The West Texas Mesonet 
was initiated by Texas Tech University in 1999 and currently monitors soil moisture at 53 
stations at depths of 5 cm, 20 cm, 60 cm, and 75 cm using water content reflectometers (615, 
Campbell Scientific, Inc., Logan, UT) (Schroeder et al., 2005). In addition the network monitors 
wind information, atmospheric pressure, solar radiation, soil temperature, precipitation, and leaf 
wetness. The Georgia Automated Environmental Monitoring Network began in 1991 
(Hoogenboom, 1993) and has since grown to include 81 stations. Soil moisture sensors have 
been added to these stations at a depth of 30 cm for the purpose of agricultural and 
meteorological monitoring. The newest large-scale soil moisture networks in the US are the 
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735 COSMOS and GPS-IR networks described in preceding sections of this manuscript. Additional 

736 networks are on the horizon as well, including the National Ecological Observatory Network 

737 (NEON) which will operate study sites in 20 eco-climatic domains throughout the U.S. in the 

738 coming years (Keller et al., 2008). 

739 Large-Scale Soil Moisture Networks Outside the United States 

740 In recent years, several large-scale soil moisture monitoring networks have been 

741 established outside of the US, either serving research purposes, supporting natural hazard 

742 forecasting, or being an integrative part of meteorological observing systems (e.g. Calvet et al., 

743 2007). Table 1 gives an overview of known large-scale networks that are currently measuring 

744 soil moisture on an operational or quasi-operational basis. No active network outside the US has 

745 a spatial extent as large as that of the US national networks, but several have spatial extents and 

746 densities comparable to the state level networks in the US. Worth mentioning are the networks, 

747 such as those in France and Mongolia, that were installed for validating satellite soil moisture 

748 missions, and thus have a setup that allows for representing as accurately as possible soil 

749 moisture variations at the spatial scale of a satellite footprint. 

750 The networks described in the previous section have each been designed to meet different 

751 research and operational objectives, and this has resulted in a large variety of measurement 

752 setups and techniques, available metadata, data access points, and distribution policies. The first 

753 action to offer a centralized access point for multiple, globally available in-situ soil moisture data 

754 sets was the Global Soil Moisture Data Bank (Robock et al., 2000; Robock et al., 2005). The 

755 GSMDB collected data sets existing at that time but did not perform any harmonization of 

756 variables or data formats. The first international initiative addressing the latter has been 

757 FLUXNET (Baldocchi et al., 2001), a “network of networks” dedicated to monitor land- 
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atmosphere exchange of carbon, energy, and water. Unfortunately, within FLUXNET soil 
moisture is not measured at all sites while, more importantly, practical use of soil moisture data 
from FLUXNET is severely hampered by restricted accessibility and the large time gap between 
acquisition of the data and making them available to the science community. 

In 2009, the International Soil Moisture Network (ISMN; http://ismn.geo.tuwien.ac.at/) 
was initiated to overcome the issues of timeliness in data delivery, accessibility, and 
heterogeneity of data (Dorigo et al., 2011a; Dorigo et al., 2011b). This international initiative is a 
result of the coordinated efforts of the Global Energy and Water Cycle Experiment (GEWEX) in 
cooperation with the Group of Earth Observation (GEO) and the Committee on Earth 
Observation Satellites (CEOS) to support calibration and validation of soil moisture products 
from remote sensing and land surface models, and to advance studies on the behavior of soil 
moisture over space and time. The decisive financial incentive was given by the European Space 
Agency (ESA) who considered the establishment of the ISMN critical for optimizing the soil 
moisture products from the SMOS mission. 

The ISMN collects and harmonizes ground-based soil moisture data sets from a large 
variety of individually operating networks and makes them available through a centralized data 
portal. Currently, the database contains almost 6000 soil moisture data sets from almost 1500 
sites, distributed among 37 networks worldwide (Fig. 14). Also the data sets contained in the 
former GSMDB were harmonized and transferred into the ISMN. It should be noted that not all 
networks are still active. 

Recently, several updates of the ISMN system were performed to keep up with the 
increasing data amount and traffic, and to meet the requirements of advanced users. Many 
datasets from operational networks (e.g., SCAN, the US Climate Reference Network, SWEX 
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781 Poland, and ARM) are now assimilated and processed in the ISMN on a fully automated basis in 

782 near-real time. In addition, an enhanced quality control system is currently being implemented 

783 (Dorigo et al., 2013) while novel methods are being explored to obtain objective measures of 

784 reliability and spatial representativeness of the various sites (Gruber et al., 2013). 

785 The steadily increasing number of soil moisture monitoring stations goes hand in hand 

786 with the growing awareness of the role of soil moisture in the climate system. Nevertheless, Figs. 

787 14 and 15 show that the current stations are concentrated geographically and mainly represent a 

788 limited number of climate classes in temperate regions. The number of permanent soil moisture 

789 stations is still very limited in the tropics (A category), dry areas (Bw classes), and in high 

790 latitude areas (Dfc and E classes). Especially in the latter the hydrological cycle is not yet well 

791 understood, and these regions are expected to be particularly sensitive to climate change. Thus, 

792 international efforts should concentrate on expanding networks in these areas. 

793 However, the major challenge is not only to setup new networks but also to keep them 

794 operational in the future. Since many networks heavily rely on project funding, their continuation 

795 is typically only guaranteed for the lifetime of the project. Thus, internationally coordinated 

796 effort should focus on developing mechanisms for continued financial and logistical support. 

797 One of such mechanisms may be the development of a soil moisture component as part of the 

798 global terrestrial network for hydrology (GTN-H) envisaged by the GCOS (GCOS, 2010). The 

799 task of such a network should go beyond the achievements of the ISMN and also define 

800 standards for the measurements themselves in order to guarantee the consistency between sites. 

801 Also, the integration of soil moisture monitoring sensors in existing operational meteorological 

802 stations would increase the probability for continued operation. 

803 
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APPLICATIONS OF LARGE-SCALE SOIL MOISTURE MEASUREMENTS 

Drought Monitoring 

Droughts are typically classified as either meteorological, agricultural, or hydrological 
(Mishra and Singh, 2010). Meteorological drought is indicated by a lack of precipitation over a 
specified region during a particular period of time. Agricultural drought occurs when declining 
soil moisture levels negatively impact agricultural production. Some have used the term 
“ecological drought” to designated similar conditions which reduce primary productivity in 
natural ecosystems (Le Houerou, 1996). These two drought concepts are closely related and 
should perhaps be represented by the composite term “agro-ecological drought”. A third 
common drought classification is hydrological drought which is a period of inadequate surface 
and subsurface water resources to support established water uses. Soil moisture is most directly 
related to agro-ecological drought which is often preceded by meteorological drought and comes 
before hydrological drought. This places soil moisture squarely in the center of the spectrum of 
drought classifications and drought indicators, but soil moisture measurements have been largely 
neglected in the science and practice of drought monitoring to date. 

In earlier decades this deficiency was unavoidable because sufficient soil moisture data 
were simply not available to enable their use in operational drought monitoring. That situation 
began to change dramatically in the 1990s with the rise of large-scale soil moisture monitoring 
networks in the US (Hollinger and Isard, 1994; McPherson et ah, 2007; Schaefer et ah, 2007), a 
change now spreading around the world. Even more recently, global maps of surface soil 
moisture based on satellite remote sensing have become available, and these could be useful in 
drought monitoring. The primary impediment to the use of soil moisture measurements in 
operational drought monitoring is no longer a lack of data, but rather a lack of scientific 
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827 understanding regarding how soil moisture measurements quantitatively indicate agro-ecological 

828 drought. Strong and transparent conceptual models are needed to link soil moisture 

829 measurements with vegetation impacts in agricultural and ecological systems. 

830 The first known attempt to use large-scale soil moisture measurements in drought 

831 monitoring was the Soil Moisture Index (SMI) introduced by Sridhar et al. (2008) based on data 

832 from the Automated Weather Data Network (AWDN) in Nebraska. Their results demonstrated 

833 that soil moisture data from 37 stations in Nebraska provided a strong quantitative drought 

834 indicator. The SMI was subsequently revised by Hunt et al. (2009) who proposed the following 

835 relationship 

836 SMI = -5 + \0F aw [2] 

837 where F A w is the fraction of available water. Fraction of available water is calculated by 

838 F AW = (0 - 9fc )/(6fc - 9wp ) [3] 

839 where 0 is the volumetric water content at a specified depth, 0f c is the volumetric water content 

840 corresponding to field capacity, and 0wp is the volumetric water content corresponding to 

841 permanent wilting point. The use of F A w as the basis for SMI is substantiated by current 

842 scientific understanding of plant water stress because water stress is more strongly related to the 

843 relative amount of plant available water in the soil than to the absolute amount of soil moisture 

844 (Allen et al., 1998). Values of F AW are typically between 0 and 1, however both higher and 

845 lower values are possible. The scaling relationship in Eq. [2] thus causes SMI values to typically 

846 fall in the range from -5 to +5. This scaling was perhaps chosen to make the range of SMI 

847 comparable to the range of other drought indicators (e.g. Palmer Drought Severity Index; Palmer, 

848 1965). Although stress thresholds vary somewhat with plant species and weather conditions, 

849 generally F AW values <0.5 result in water stress (Allen et al., 1998). When F AW is 0.5, the SMI 
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value is 0, the transition between stressed and non-stressed conditions. Again using data from 
the Nebraska AWDN, Hunt et al. (2009) found that the modified SMI was effective for 
identifying drought onset as well as soil recharge from rainfall events following significant dry 
periods. Recently, the SMI was applied using data from a network of six monitoring stations in 
the Czech Republic (Mozny et al., 2012). That study supported the drought intensity scheme 
proposed by Sridhar et al. (Sridhar et al., 2008) in which SMI values <-3 signify severe or 
extreme drought. Mozny et al. (Mozny et al., 2012) related the concept of “flash drought” to 
SMI, specifying that a flash drought occurs when SMI values decrease by at least 5 units during 
a period of 3 weeks. Thus, the SMI concept has shown good potential as a quantitative drought 
indicator based on soil moisture measurements, but some key uncertainties remain. The 
indicator is sensitive to the site- and depth-specific values chosen for 0f c and 0wp. These critical 
water contents can be estimated from the in situ soil moisture time series in some cases (Hunt et 
al., 2009), measured directly in the laboratory, calculated using pedotransfer function models 
(Schaap et al., 2001), or estimated from literature values (Sridhar et al., 2008), but best practices 
for determining these parameters in the SMI context need to be developed. 

Recently, Torres et al. (2013) introduced a method for using long-term measurements of 
soil water deficits (SWD) from a large-scale monitoring network to compute site-specific 
drought probabilities as a function of day of year. Improved quantification of seasonal patterns 
in drought probability may allow crop cycles to be better matched with periods when drought is 
less likely to occur; therefore, yield losses due to drought may be reduced. Soil water deficit for 
each soil layer (D) is defined as 

D = (e fc -0)Az [4] 
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872 where Az is the thickness of the soil layer, and SWD is calculated by summing D over the 

873 desired soil layers. Soil moisture data from eight stations of the Oklahoma Mesonet spanning 15 

874 years were used to calculate deficits for the 0-10 cm, 10-40 cm, and 40-80 cm layers. Drought 

875 was defined in this context as a period when SWD is sufficient to cause plant water stress, i.e. 

876 SWD exceeds a predetermined threshold. The threshold was set for each site and layer as 

877 0.5TAW, where TAW is the total available water calculated by substituting 0 wp for 0 in Eq. [4], 

878 Values of SWD calculated from 0-40 cm (SWD 40 ) were similar to 7-d cumulative atmospheric 

879 water deficits (AWD), calculated as reference evapotranspiration minus precipitation, during 

880 much of the spring and fall, but the soil and atmospheric deficits diverged in the winter and 

881 summer months (Fig. 16). 

882 Historical drought probabilities estimated for each day of the year using the SWD data 

883 were consistent between depths and agreed with general knowledge about the climate of the 

884 region (Fig. 17), while probabilities estimated using AWD data (Purcell et al., 2003) were 

885 substantially lower and inconsistent with general knowledge about the region and with prior 

886 drought probability estimates in nearby states. Torres et al. (2013) proposed modifications to the 

887 AWD method, either lowering the AWD threshold used to define drought or extending the 

888 summation period from 7 to 15 days, both of which resulted in drought probability estimates 

889 more consistent with the estimates from SWD method. They concluded that the new SWD 

890 method gave plausible and consistent results when applied to both the 0- to 40- and 0- to 80-cm 

891 soil layers and should be utilized when long-term soil moisture data are available. 

892 The first known operational use of large-scale soil moisture measurements for drought 

893 monitoring involves, not SMI or SWD, but a related measure, plant available water (PAW). 

894 Plant available water is defined as 
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^ (3 "3m)*, [51 

for soil layers i= 1 ...n of thickness dz,. In 2012, the Oklahoma Mesonet (McPherson et al., 2007) 
introduced daily-updated PAW maps based on its network of >100 stations monitoring soil 
moisture at standard depths of 5, 25, and 60 cm. These maps are intended for use in drought 
monitoring and show PAW for the 0-10 cm (4-inch), 0-40 cm (16-inch), and 0-80 cm (32-inch) 
soil layers (http://www.mesonet.org/index.php/weather/category/soil_moisture_temperature). 

The depth units (e.g. mm or inches) of PAW make it compatible with familiar hydrologic 
measurements such as precipitation and evapotranspiration. Figure 18 shows maps of total 
rainfall, total short-crop reference ET based on the FAO-56 procedure (Allen et al., 1998), and 
average PAW across the state of Oklahoma during May 2012. Dry conditions prevailed across 
the state with reference ET exceeding rainfall at all measured locations. The PAW map reflects 
the influence of rainfall and ET with relatively high PAW values in eastern, northeastern, and 
central OK corresponding to regions with relatively high rainfall and/or relatively low reference 
ET. However, the PAW maps also suggest more complex influences of vegetation, soil type, 
and landscape “memory”. For example, note that PAW values were generally lower in the 
southwest portion of the state than in the Panhandle region even though the Panhandle region 
experienced lower rainfall totals and comparable reference ET. This illustrates the challenges 
with using atmospheric data alone to monitor agro-ecological drought and suggests a unique and 
complementary role for soil moisture measurements. 

These recent developments in the use of soil moisture measurements for drought 
monitoring are encouraging; however the research needs in this area are significant. As yet, little 
is known regarding how soil moisture-based drought indicators relate to other widely-accepted 
drought indicators like the Standardized Precipitation Index (Guttman, 1999) or the Palmer 


40 

5585 Guilford Rd„ Madison Wl 53711 


Page 41 of 110 


Soil Science Society of America Journal 


918 Drought Severity Index. Likewise, we do not know how soil moisture-based drought indicators 

919 are related to actual drought impacts in agricultural or ecological systems. Already SMI, SWD, 

920 and PAW have demonstrated potential as soil moisture-based drought indicators driven by in situ 

921 measurements, but these three indicators all address the question, “How dry is it?” rather than the 

922 equally important question, “How much drier than average is it?” Other soil moisture-based 

923 indicators have been proposed on the basis of numerical modeling studies. These include the 

924 model-based Normalize Soil Moisture index (Peled et al., 2010) and the Soil Moisture Deficit 

925 Index (Narasimhan and Srinivasan, 2005), neither of which has been evaluated using actual soil 

926 moisture measurements. 

927 Furthermore, most in situ soil moisture measurements are made under grassland 

928 vegetation because of problems with establishing long-term meteorological stations in cropland 

929 or forest. There is a dearth of research on how to estimate soil moisture under contrasting land 

930 use/land cover combinations based on in situ observations under grassland vegetation. This 

931 deficiency complicates the interpretation of agro-ecological drought indicators based on in situ 

932 soil moisture measurements. Clearly, there should be a role for satellite remote sensing of soil 

933 moisture to assist in overcoming some of the deficiencies of drought monitoring by in situ soil 

934 moisture observations. Bolten et al. (2010) showed that AMSR-E surface soil moisture retrievals 

935 could add significant value to root zone soil moisture predictions in an operational drought 

936 modeling framework. However, there are as yet no operational systems for drought monitoring 

937 that utilize satellite soil moisture measurements. We anticipate a surge in this type of research in 

938 the near future. 

939 

940 Meteorological Modeling and Forecasting 
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Drought provides a clear example of the interaction between the atmosphere and the land 
surface, an interaction strongly influenced by the soil moisture conditions. A schematic of 
atmospheric boundary layer (ABL) interactions with the land surface is presented in Fig. 19. 
Daytime growth of the ABL is directly affected by soil and vegetation states and processes, and 
these processes play a role in partitioning the energy balance which relates net radiation to soil 
heat flux, sensible heat flux, and latent heat flux, i.e. evapotranspiration. Root-zone soil moisture 
can influence the atmospheric boundary layer by controlling land surface energy and moisture 
fluxes. For example, Basara and Crawford (2002) found that soil water content in the root-zone 
during the summer was linearly correlated with daytime evaporative fraction and daily- 
maximum values of sensible heat flux and latent heat flux on days with strong radiative forcing 
and weak shear in the lower troposphere. Root-zone soil moisture was also linearly related to 
key parameters in the ABL such as the daily maximum air temperature at 1.5 m. 

Numerous large-scale hydrologic-atmospheric-remote sensing experiments have been 
conducted to better understand the soil moisture-moderated interactions of the soil-vegetation 
system with the diurnal atmospheric boundary layer. Improved parameterization of general 
circulation models (GCMs) was one of the initial objectives of the experiments. Table 2 gives a 
concise overview of a few of these experiments, including HAPEX-MOBILHY which was the 
first experiment conducted on this scale (Andre et al, 1986; Andre et al., 1988). It should be 
noted most of the experiments listed cover large geographic areas which play a significant role in 
the general circulation system of the planet. 

The strong linkage of surface soil moisture and parameterization of soil hydraulic 
processes with ABL response was demonstrated by Ek and Cuenca (1994), based on data from 
the HAPEX-MOBILHY. This study found that variations in soil hydraulic process 
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964 parameterization could have a clear impact on the simulated surface energy budget and 

965 atmospheric boundary layer development. This impact was accentuated for dry to moderate soil 

966 moisture conditions with bare soils. Ek continued to do considerable work in the area of 

967 simulation of the ABL and the influence of soil moisture conditions, often using data from 

968 regional experiments such as HAPEX-MOBILHY and the Cabauw data set from the Netherlands 

969 (Monna and van der Yliet, 1987). Data from HAPEX-MOBILHY was used to evaluate the 

970 evolution of the relative humidity profile in the ABL in Ek and Mahrt (1994). The relationships 

971 between canopy conductance, root density, soil moisture and soil heat flux with simulation of the 

972 ABL using the Cabauw data set was investigated in Ek and Holtslag (2004). It should be noted 

973 that the ABL simulation evolved from the Oregon State University 1-D planetary boundary-layer 

974 model (OSU1DPBL) (Mahrt and Pan, 1984; Pan and Mahrt, 1987) to the Coupled Atmospheric 

975 boundary layer-Plant-Soil (CAPS) model. These models in turn are the basis for the Noah land- 

976 surface model (Chen and Dudhia, 2001; Ek et ah, 2003) which plays a major role in the Medium- 

977 Range Forecast Model for numerical weather prediction (NWP) at the NOAA National Center 

978 for Environmental Prediction. 

979 Given its influence on ABL development, root-zone soil moisture can have a strong 

980 influence on weather forecasts. If not suitably constrained, the root-zone soil moisture in a NWP 

981 model will drift from the true climate, resulting in erroneous boundary layer forecasts (Drusch 

982 and Viterbo, 2007). Root-zone soil moisture cannot currently be observed at the spatial scales 

983 required by NWP, and since the mid 1990s many NWP centers have been indirectly constraining 

984 their model soil moisture using methods that minimize the errors in measured screen-level (1.5- 

985 2.0m) temperature and humidity (Best et al 2007; Hess, 2001; Mahfouf 1991; Mahfouf et al 

986 2009). While this approach reduces boundary layer forecast errors, it does not generate realistic 
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soil moisture since the latter is often adjusted to compensate for model errors unrelated to soil 
moisture (Douville et al., 2000; Drusch and Viterbo, 2007; Hess, 2001). Ultimately a model with 
inaccurate soil moisture cannot accurately describe the atmosphere across the full range of 
forecast lengths produced from NWP models. 

Hence, the NWP community has been working towards improving model soil moisture 
by assimilating remotely sensed near-surface soil moisture. Near-surface soil moisture is more 
directly related to root-zone soil moisture than screen-level variables, and assimilating near- 
surface soil moisture data has been shown to improve model root-zone soil moisture (Calvet et 
al., 1998; Hoeben and Troch, 2000; Montaldo et al., 2001). Figure 20 compares several 
experiments constraining model root-zone soil moisture by assimilating observations of near- 
surface soil moisture and screen-level temperature and relative humidity, highlighting the 
fundamental difference between these two approaches. These experiments were conducted with 
Meteo-France’s NWP land surface model using an Extended Kalman Filter and the Advanced 
Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) Land Parameter 
Retrieval Model near-surface soil moisture data (Owe et al., 2008). Refer to Draper et al. (2011) 
for further details. 

In general, assimilating the screen-level observations (black-dashed line) improved the fit 
between the mean forecast and observed screen-level variables, compared to the open loop (no 
assimilation, solid black line). However, the assimilation had a slight negative impact on the fit 
between the mean forecast and observed near-surface soil moisture. In contrast, assimilating the 
AMSR-E soil moisture (grey solid line) improved the fit between the mean forecast and observed 
near-surface soil moisture, while degrading the fit between the modeled and observed screen- 
level variables. This result is consistent with previous studies showing that adjusting model soil 


44 

5585 Guilford Rd., Madison Wl 53711 


Page 45 of 110 


Soil Science Society of America Journal 


1010 moisture to improve screen-level forecasts does not necessarily improve soil moisture (Douville 

1011 et al, 2000; Drusch and Viterbo, 2007; Seuffert et al 2004), and conversely improving the model 

1012 soil moisture does not necessarily improve atmospheric forecasts (Seuffert et al 2004). 

1013 Consequently, in the foreseeable future it is unlikely that remotely sensed near-surface soil 

1014 moisture will be used in NWP in place of screen-level observations. However, combining the 

1015 assimilation of both observation types can reduce errors in both model soil moisture and low- 

1016 level atmospheric forecasts. For example, when both data types were assimilated together (grey 

1017 dashed line) in Fig. 20, the fit between the model and both observation types was improved, 

1018 although the mean soil moisture improvements were very small (see also Seuffert et al, 2004). 

1019 Currently near-surface soil moisture observations are assimilated operationally at the UK 

1020 Met Office (UKMO) and the European Centre for Medium Range Weather Forecasting 

1021 (ECMWF). While the development of soil moisture assimilation in NWP is motivated by the 

1022 eventual use of L-band observations (e.g. SMOS and SMAP), both centers are assimilating 

1023 Advanced Scatterometer (ASCAT) Surface Degree of Saturation (SDS) data (Bartalis et al, 

1024 2007), since this is currently the only operationally-supported remotely sensed soil moisture 

1025 product. At the UKMO the screen-level observation based soil moisture analysis was amended 

1026 in July 2010 to also constrain the near-surface soil moisture by nudging it with ASCAT SDS data 

1027 (Dharssi et al, 2011). Compared to nudging with only screen-level observations, adding the 

1028 ASCAT data very slightly improved near-surface soil moisture forecasts over selected sites in 

1029 the US, while also improving screen level temperature and relative humidity forecasts over the 

1030 tropics and Australia (with neutral impact elsewhere). At ECMWF the NWP land surface 

1031 analysis was updated in November 2010 to an extended Kalman filter based scheme, enabling 

1032 the assimilation of remotely sensed data (de Rosnay et al 2012, Drusch et al 2009). The ASCAT 
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SDS are not used in their weather forecasting model, but are assimilated together with screen- 
level observations in an offline land surface analysis system. Including the ASCAT data in this 
system has had a neutral impact on near-surface soil moisture and screen-level forecasts 
(Albergel et al 2012b; de Rosnay et al, 2012). 

The above examples highlight some challenges of land data assimilation specific to NWP 
applications. For example, the computation cost of the assimilation is a major limitation in NWP 
(de Rosnay et al 2012, Drusch et al 2009), hence the assimilation methods applied must be 
relatively simple. Further work is required to improve the land surface analysis schemes used in 
NWP, and in particular to propagate the surface soil moisture information into the root-zone (not 
currently achieved by the schemes in place at the UKMO or ECMWF). Additionally, Dharssi et 
al. (201 1) and de Rosnay et al. (2012) identified the observation bias correction strategy, i.e. the 
method by which satellite derived surface soil moisture values are adjusted to be consistent with 
the model used for assimilation, as a likely cause of the limited impact of assimilating the 
ASCAT data. Bias correction of remotely sensed soil moisture is difficult in NWP, since the long 
data records required to estimate statistics of the model climatology are not available from NWP 
models, due to frequent model updates and the prohibitive cost of rerunning models. 

However, the greatest challenge faced by soil moisture assimilation in NWP is that 
improving the model soil moisture may not immediately improve atmospheric forecasts, due to 
errors in the model physics. It is likely that the greatest contribution of using remotely sensed 
near-surface soil moisture observations in NWP will be in helping to identify and address these 
physics errors. Already, the availability of remotely sensed soil moisture and efforts to 
assimilate that data have stimulated improvements in modeling soil moisture processes. For 
example, in response to discrepancies between modeled and SMOS observed brightness 
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1056 temperatures, ECMWF recently improved their bare soil evaporation parameterization, resulting 

1057 in improved model near-surface soil moisture and brightness temperature (Albergel et al, 2012b). 

1058 As soil moisture data is used more extensively in NWP models, this should also help to expose 

1059 and eventually address other errors in the model surface flux processes. 

1060 Ecological Modeling and Forecasting 

1061 Ecological modeling is another area which could logically benefit from increased 

1062 availability of large-scale soil moisture monitoring. Soil moisture is a key parameter in the 

1063 control of plant growth, soil respiration, and distribution of plant functional types in terrestrial 

1064 ecosystems (Blyth et al. 2010; Ren et al. 2008; Pan et al. 1998; Neilson 1995). Plant growth (i.e. 

1065 assimilation of C02 through photosynthesis) is coupled with water loss through transpiration 

1066 which is regulated by soil water availability (Yang et al. 2011; Sellers et al. 1997; Field et al., 

1067 1995). Decomposition of soil organic carbon is also sensitive to soil moisture content via 

1068 microbial activity and other processes (Ise and Moorcroft 2006; Xu et al. 2004; Orchard and 

1069 Cook 1983). Furthermore, temporal and spatial availability of soil moisture content constrains 

1070 distribution and properties of plant functional types (Bremond, Boom, and Favier 2012; 

1071 Seneviratne et al. 2010; Gerten et al. 2004; Breshears and Barnes 1999). 

1072 A striking example of the interactions between vegetation and soil moisture conditions is 

1073 provided by the Tiger Bush sites in the HAPEX-Sahel experiment. The Tiger Bush is made up 

1074 of relatively long and narrow patches of vegetation approximately 40-m wide separated by 

1075 nearly cemented patches of bare soil approximately 60-m wide and these sites are characteristic 

1076 of certain regions in the Sahel. One can note in the >3-m deep profile in Fig. 21 (monitored by 

1077 neutron probe) that there is limited variation in soil moisture content and only in the upper 50 cm 

1078 of the bare soil profile, while there are appreciable soil moisture changes even past 300-cm in the 
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vegetated strip. The result is that nearly all of the high intensity rainfall during the rainy season 
in this environment runs off the bare soil into the vegetated strip which thereby receives on the 
order of two hundred percent of the precipitation. Yerhoef (1995) noted this effect and that the 
result was a well- watered vegetation strip adjacent to a very dry bare soil strip in this 
environment. Verhoef (1995) was able to show that in the generally hot and dry conditions of 
the Sahel, advective conditions for sensible heat flux from the bare soil resulted such that the 
evapotranspiration (ET) from the vegetated strip clearly exceeded the potential, or reference, ET 
rate (Verhoef et al., 1999; Verhoef and Allen, 2000). Carbon fluxes would obviously be affected 
by the heterogeneity in the Tiger Bush system, as well. 

To better understand and predict ecosystem dynamics such as these, different classes of 
ecological models have been developed for various scales and emphases. For example, 
biogeography models such as MAPSS (Neilson, 1995) and BIOMEs (Prentice et al., 1992; 
Haxeltine and Prentice, 1996) focus on the distribution of species and ecosystems in space. 
Biogeochemistry models such as CENTURY/DAYCENT (Parton et al., 1987, 1998), RothC 
(Jenkinson and Coleman, 1994) and DNDC (Li et al., 1992) place emphasis on the carbon and 
nutrient cycles within ecosystems. Biophysics models based on soil-vegetation-atmosphere 
transfer (SVAT) schemes (SiB: Sellers et al. 1986; BATS: Dickinson et al. 1986) have been 
developed to support regional and global climate modeling to provide accurate information for 
the fluxes of water, radiation, heat and momentum between the atmosphere and the various land 
surfaces. Recently developed dynamic global vegetation models (DGVM) such as LPJ (Sitch et 
al., 2003), IBIS (Foley et al., 1996) and MCI (Bachelet et al., 2001), generally incorporate above 
classes of models and schemes to simulate dynamics of potential vegetation and its associated 
biogeochemical and hydrological cycles. 
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1102 These models estimate soil moisture content or its proxy using different schemes such as 

1103 the bucket method (Robock et al. 1995; Manabe 1969), the precipitation to potential 

1104 evapotranspiration ratio method (Scheffer et al., 2005), and the water balance model (Law et al. 

1105 2002). Details of these and other schemes are discussed in Shao and Henderson-Sellers (1996) 

1106 and Ren et al. (2008). These schemes often use simple algorithms to reduce computational 

1107 demand and are less reliable compared to schemes used in hydrologic models [e.g. the Richards 

1108 equation (Richards, 1931)]. Also, especially in cases of large scale ecological models, a more 

1109 realistic parameterization of soil moisture content at subgrid-scale as related to topography is 

1110 often desirable (Gordon et al. 2004). Optimization of the degree of the simplification and the 

1111 spatial resolution is necessary due to computational restrictions, but is difficult to judge due to 

1112 lack of adequate observational data with which to verify the performance of the models (Ren et 

1113 al.2008). 

1114 Traditionally, ecological models have been tested through intercomparison studies such 

1115 as the Vegetation/Ecosystem Modeling and Analysis Project (YEMAP; YEMAP Members 

1116 1995), the Carbon Land Model Project (CLAMP; Randerson et al. 2009), the Project for 

1117 Intercomparison of Land-surface Parameterization Schemes (PILPS; Henderson-Sellers et al. 

1118 1996; 1995), and the Global Soil Wetness Projects (GSWP/GSWP2; Dirmeyer et al. 2006; 

1119 Dirmeyer 1999) because evaluating the model performance, especially at larger scales, has been 

1120 difficult due to the incompleteness of observation data sets. However, these models are not 

1121 independent because they have integrated the same theories and relied on similar data sets as 

1122 they evolved (Reichstein et al. 2003). Therefore, while model intercomparison is an important 

1123 task, extreme care must be exercised to derive any conclusions. 
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Future research advances in this area will require use of newly available observation data 
(Seneviratne et al. 2010). Observation data from large-scale soil moisture monitoring in 
particular should be valuable to validate the simplification and scaling of ecological models. 
Wagner et al. (2003) found that modeled soil moisture from the Lund-Potsdam-Jena (LPJ) 
dynamic global vegetation model agreed “reasonably well” over tropical and temperate locations 
with values derived from satellite microwave scatterometer, yielding Pearson correlation 
coefficients >0.6. The agreement was poorer over drier and colder climatic regions. However, 
few studies have used large-scale soil moisture data to improve the structure or parameterizations 
of ecological models or to improve model predictions through data assimilation. 

Furthermore, the relationship between soil moisture and terrestrial ecosystem is dynamic 
and interdependent: soil moisture constrains the properties of the ecosystem as described earlier, 
which in turn, modifies hydrology through evapotranspiration, LAI, and surface roughness 
(Breshears and Barnes 1999). Newer generations of ecological models, especially dynamic 
global vegetation models, include these important feedback processes to simulate the effects of 
future climate change on natural vegetation and associated carbon and hydrologic cycles. 
Validation of these models studies will reveal an increased need for data from large-scale soil 
moisture observations across various ecosystems and for continuous expansion of observation 
networks. 

Hydrologic Modeling and Forecasting 

One motivation underlying many large-scale soil moisture monitoring efforts is the desire 
to more accurately model and forecast watershed dynamics, especially streamflow and flood 
events. Pauwels et al. (2001) demonstrated the possibility of improved stream discharge 
estimates through assimilation of surface soil moisture estimates derived using data from the 
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1147 ESA satellites ERS1 and ERS2 into a land atmosphere transfer scheme. The study was limited 

1148 to bare soil conditions and small catchments (<20 km ). The assimilation improved discharge 

1149 estimates 20-50% in seven out of 12 cases considered, but degraded model accuracy by up to 

1150 10% in the remaining five cases. Francois et al. (2003) showed that assimilation of ERS1 SAR 

1151 data into a simple two-layer land surface scheme through an extended Kalman filter approach 

1152 improved the Nash-Sutcliffe efficiency (NSE) for streamflow from 70% to 85%. Their study 

1153 involved a larger catchment (104 km ) than that of Pauwels et al. (2001) and included vegetation 

1154 cover. The sensitivity of simulated flow to soil moisture was highest when soil moisture itself 

1155 was high. The assimilation scheme was also able to correct for 5-10% errors in the input data, 

1156 e.g. potential evapotranspiration or precipitation. 

1157 More recent applications of large-scale soil moisture data for hydrologic modeling and 

1158 forecasting have focused on newer satellite remote sensing datasets. Brocca et al. (2010) used a 

1159 simple nudging scheme to assimilate the ASCAT surface soil wetness index into a rainfall — 

1160 runoff model for five catchments (<700 km ) in the Upper Tiber River basin in Italy. 

1161 Assimilation increased the NSE for streamflow prediction during flood events in all five 

1162 catchments, with increases ranging from 2 to 17% (Fig. 22). In a subsequent study, root zone 

1163 soil moisture was estimated from the ASCAT surface soil moisture data through application of 

1164 an exponential filter, and both data types were then assimilated into a two-layer rainfall — runoff 

1165 model using an ensemble Kalman filter approach (Brocca et al., 2012). Assimilation of the root 

2 

1166 zone soil moisture estimates produced a clear improvement in discharge prediction for a 137 km 

1167 catchment (NSE improved from 76% to 86%), while assimilation of surface soil moisture had 

1168 only a small effect. 


51 

5585 Guilford Rd„ Madison Wl 53711 


1169 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 

1181 

1182 

1183 

1184 

1185 

1186 

1187 

1188 

1189 

1190 


Soil Science Society of America Journal 


Page 52 of 110 


Thus far only a few studies have evaluated methods for using soil moisture data to 
improve hydrologic forecasting in catchments of >1000 km . One example is the work of Meier 
et al. (2011) in which the ERS1 and ERS2 soil water index was used, along with rainfall data, to 
drive a conceptual rainfall — runoff model in an ensemble Kalman filter framework assimilating 
observed discharge every 10 days. The method was applied to three catchments in the Zambezi 
River basin in southern Africa. The catchments ranged in size from 95,300 to 281,000 km . The 
catchment average soil water index correlated well with measured discharge when the data were 
shifted by a catchment-specific time lag. This time lag allowed 40-d lead time streamflow 
forecasts with a NSE value of 85% for the largest watershed, but in a catchment with steep 
slopes and little soil water storage the lead time was as short as 10 d. Gains in streamflow 
forecast accuracy, especially during flood events, have even been demonstrated by using point 
soil moisture observations from a single location within a catchment of 1 120 km , suggesting 
that even sparse observations networks in large catchments can be quite useful (Aubert et al., 
2003). 

That large-scale soil moisture monitoring can benefit hydrologic modeling and 
forecasting is now well-established with gains in forecast efficiency of 10-20% being typical; 
however, significant challenges and uncertainties remain. Most of the research to date in this 
area has focused on the use of satellite derived surface soil moisture products, with few studies 
using in situ soil moisture measurements within a data assimilation framework (Aubert et al., 
2003; Chen et al., 2011). Thus, the world’s growing in situ soil moisture monitoring 
infrastructure (Table 1) is a virtually unexplored resource in this context, and many opportunities 
exist to develop hydrologic forecasting tools which utilize that infrastructure. 
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1191 A key challenge associated with assimilation of soil moisture data, regardless of the 

1192 source, is to identify and overcome structural deficiencies in the hydrologic models themselves. 

1193 For example, a data assimilation experiment using in situ soil moisture measurements in 

1194 Oklahoma was generally unsuccessful in improving streamflow predictions from the widely used 

1195 Soil and Water Assessment Tool (SWAT) model (Chen et al., 2011). The calibrated SWAT 

1196 model significantly underestimated the vertical coupling of soil moisture between upper and 

1197 lower soil layers, and the inadequate coupling was apparently a structural, rather than parametric, 

1198 problem in the model. Thus, the ensemble Kalman filter assimilation approach was not effective 

1199 in improving estimates of deep soil moisture or streamflow. This particular challenge of correctly 

1200 representing linkages between soil moisture across two or more soil layers has been identified as 

1201 a key concern in studies with other models as well (Brocca et al., 2012). Further research is 

1202 needed to optimize the structure of SWAT and other hydrologic models in order to maximize the 

1203 benefits from assimilating increasingly available large-scale soil moisture observations (Brocca 

1204 etal.,2012). 

1205 Another challenge which has been encountered in this arena is uncertainty regarding 

1206 proper characterization of model errors and observation errors within the assimilation procedure 

1207 (Francois et al., 2003; Brocca et al., 2012). Statistical representations of model errors must often 

1208 be made in a somewhat arbitrary or subjective fashion, and pre-existing biases in either the 

1209 observations or the model can be particularly problematic (Chen et al., 2011; Brocca et al., 

1210 2012). Nevertheless, research in this area appears to be gaining momentum, and opportunities 

1211 abound to advance hydrologic modeling and forecasting with the help of existing and emerging 

1212 large-scale soil moisture datasets. 

1213 
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PRIMARY CHALLENGES AND OPPORTUNITIES 

In this review, we have attempted to describe the state of the art in large-scale soil 
moisture monitoring and to identify some critical needs for research to optimize the use of 
increasingly available soil moisture data. We have considered: 1) emerging in situ and proximal 
sensing techniques, 2) dedicated soil moisture remote sensing missions, 3) soil moisture 
monitoring networks, and 4) applications of large-scale soil moisture measurements. The 
primary challenges and opportunities in these topic areas can be summarized as follows. For 
emerging in situ and proximal sensing techniques (e.g. COSMOS and GPS) empirical 
confirmations of theoretical predictions regarding the variable measurement depths are needed. 
Calibration procedures for these methods, as well as the DTS methods, need to be further refined 
and standardized with due accounting for site-specific factors such as soil and vegetation 
characteristics which influence instrument performance. Spatial and temporal heterogeneity in 
these site-specific factors must also be considered in some instances, creating additional 
challenges. Also, the community of expertise for these methods, that is the number of 
researchers with the capability to advance these technologies, needs to be expanded. 

Probably the largest share of scientific resources in this general topic area is currently 
devoted to the advancement of satellite remote sensing approaches for soil moisture monitoring. 
These investments are bearing fruit, but challenges and opportunities remain. One significant 
challenge is to further improve methods for estimating root zone soil moisture, the information 
we often need, using surface soil moisture observations, the information satellites provide. 
Progress has been made towards this goal, by using data assimilation into numerical models to 
retrieve root-zone soil moisture from near-surface observations. Continued improvements are 
also needed in calibration and validation of remotely-sensed soil moisture products because the 
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1237 relatively coarse resolution of these products is not well matched with most in situ observations. 

1238 Although not primarily a scientific challenge, more work is needed to reduce problems 

1239 associated with RFI. Similarly, continuity of missions is a necessity if remotely sensed soil 

1240 moisture data are to be adopted for operational applications like numerical weather prediction. 

1241 In contrast with remote sensing approaches, relatively few resources are currently 

1242 devoted toward large-scale in situ soil moisture monitoring networks. Although the number of 

1243 networks is growing steadily, the lack of standardization of procedures across networks is a 

1244 significant challenge. There is a need for rigorous guidelines and standards to be developed and 

1245 adopted worldwide for in situ soil moisture monitoring networks, similar to guidelines for 

1246 measurement of other meteorological variables. Best practice standards for sensor selection, 

1247 calibration, installation, validation, and maintenance are needed, as well as standards for site 

1248 selection, data quality assurance and quality control, data delivery, metadata delivery, and data 

1249 archives. The recent recognition of soil moisture as an “Essential Climate Variable” by the 

1250 Global Climate Observing System, and the development of the International Soil Moisture 

1251 Network are positive steps in this direction, but much more is needed. 

1252 For both in situ networks and remote sensing approaches, sustainability is a significant 

1253 challenge, perhaps underestimated. Societal and scientific needs for soil moisture data would 

1254 seem to justify that our monitoring systems be designed to function without interruption for 

1255 many decades. Current realities within science and society at large typically result in monitoring 

1256 systems which are designed to function for only a few years. Researchers are rewarded for 

1257 developing new systems and technologies, not for ensuring their long-term viability. Successful 

1258 long-term operation of monitoring systems generally requires substantial state or federal support. 

1259 Securing such long-term support for soil moisture monitoring systems is often difficult. Thus, 
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determining effective pathways to transition monitoring systems from research mode to 
operational mode remains a key challenge. 

In closing, we again note the growing need to develop the science necessary to make 
effective use of the rising number of large-scale soil moisture datasets. One area where 
significant progress seems possible in the near term is the use of large-scale soil moisture data 
for drought monitoring. Already some progress has been made using in situ data for this purpose 
and approaches using remote sensing data seem sure to follow. Significant efforts have been 
devoted to the use of soil moisture observations within the area of numerical weather prediction. 
In general, assimilation of soil moisture data has resulted in only modest improvements in 
forecast skill. The primary problem is that the current model structures are not well suited for 
assimilation of these data, and the model physics may not be properly parameterized to allow for 
accurate soil moisture values. A smaller effort, but arguably greater progress, has been made in 
the assimilation of soil moisture data into models designed primarily for hydrologic prediction, 
especially rainfall — runoff models. Here gains in forecast efficiency of 10-20% are not 
uncommon. Nonetheless, as with numerical weather prediction, a key challenge is to identify or 
create models that are structured in a way that is optimal for assimilation of soil moisture data. 

To date little or no progress has been made in using large-scale soil moisture observations to 
improve the structure, parameters, or forecasts of ecological models, and perhaps surprisingly, 
the same can be said for crop models. These topic areas are ripe with opportunity and challenges 
yet to be uncovered. Another frontier where the potential is great but the workers are few is the 
use of soil moisture observations in socio-economic modeling and forecasting to address issues 
such as drought impacts and food security (Simelton et al., 2012). We are optimistic that these 
challenges and opportunities can be addressed by improved communication and collaboration 
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1283 across the relevant disciplines. The international soil science community has much to contribute 

1284 in this context. Hopefully this review will be a small step towards further engaging that 

1285 community in advancing the science and practice of large-scale soil moisture monitoring for the 

1286 sake of improved Earth system monitoring, modeling, and forecasting. 

1287 
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1289 Table 1. Partial list of large-scale (>100 2 km 2 ) in situ soil moisture monitoring networks ordered 

1290 from largest to smallest in areal extent. 


Network Name 

Country or 
State 

Site 

no. 

Extent 

Density 

Reference 




km 2 

km 2 st" 1 


Inside the U.S. 






Soil Climate Analysis Network 

USA 

180 

3100 2 

230 2a 

Schaefer et al. (2007) 

Climate Reference Network 

USA 

114 

3100 2 

290 2 

Palecki and Groisman (2011) 

Cosmic Ray Soil Moisture Observing System 

USA 

67 

3100 2 

380 2 

Zreda et al. (2012) 

Plate Boundary Observatory Network 

Western US 

59 

1800 2 

240 2 

Larson et al. (2008) 

Automated Weather Data Network 

Nebraska 

53 

450 2 

62 2 

Hubbard et al. (2009) 

Oklahoma Me sonet 

Oklahoma 

108 

430 2 

41 2 

Illston et al. (2008) 

Automated Environmental Monitoring Network 

Georgia 

81 

390 2 

44 2 

Hoogenboom (1993) 

Water & Atmospheric Resources Monitoring Program 

Illinois 

19 

390 2 

89 2 

Scott et al. (2010) 

Environment and Climate Observing Network 

N. Carolina 

37 

370 2 

61 2 

Weinan et al. (2012) 

West Texas Mesonet 

Texas 

53 

300 2 

41 2 

Schroeder et al. (2005) 

ARM-SGP Extended Facilities 13 

OK/KS 

13 

150 2 

42 2 

Schneider et al. (2003) 

Outside the U.S. 






Tibet-Obs 

China 

46 

1600 2 

230 2 

Su et al. (2011) 

OzNet 

Australia 

38 

290 2 

47 2 

Smith et al. (2012) 

SMOSMANIA 

France 

21 

200 2 

44 2 

Calvet et al. (2007) 

Gourma Mesoscale Site 

Mali 

10 

170 2 

55 2 

de Rosnay et al. (2009) 

Automatic Stations for Soil Hydrology 

Mongolia 

12 

140 2 

40 2 

Yang et al. (2009) 

Central Tibetan Plateau SMTMN C 

China 

50 

100 2 

14 2 

Zhao et al. (2013) 

Umbria Region Hydro-meteorological Network 

Italy 

15 

100 2 

26 2 

www.cfumbria.it 


1291 '‘Density is calculated as the ratio of extent to site number. Note 100 2 km 2 10,000 km 2 . 

1292 b The ARM-SGP Extended Facility network is being restructured. Values listed are projections for summer 2013. 

1293 'Soil Moisture/I emperature Monitoring Network 
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1294 Table 2. Selected large-scale hydrologic-atmospheric-remote sensing experiments. 


Experiment 

Lead Agency 

Location 

Climatic Regime 

Observation Period 

HAPEX-MOBILHY 

Meteo - France 

Southwest France 

Temperate Forest 

Summer, 1986 

HAPEX-Sahel 

Meteo - France 

Niger 

Tropical Arid 

Summer, 1992 

BOREAS 

NASA 

Canada 

Boreal Forest 

Spr./Fall 1994,1996 

IHOP 

NSF 

KS, OK, TX 

Continental 

2002 

HYMeX 

GEWEX 

Europe 

Mediterranean 

2010- 2020 (LOP) 

2011- 2015 (EOP) 

CZO 

NSF 

6 sites 

Varies 

2007 - Current 

AirMOSS 

NASA 

7 sites 

Varies 

2011-2015 


1295 LOP - Long-term observation period 

1296 EOP - Enhanced observation period 

1297 


1298 


1299 

1300 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



Fig. 1 . Response function for cosmic-ray probe for soils with pore water only (solid black line) 
and those with pore water and other water, such as lattice and organic matter (dashed black line). 
N is the mesaured neutron intensity, and No is a calibration parameter representing the neutron 
intensity above dry soil. The presence of other water shifts the line horizontally from point A to 
B and A’ to B’, and the new line is steeper than the original line for the same moisture range (B- 
B’ vs. A-A’). Section B-B’ can be placed on the original line by translating it up to fall on 
section A’-A”. Thus, accounting for additional (non-pore) water does not require a new response 
function, but merely a translation along the original function by the amount equal to that non- 
pore water component. 
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21 Fig. 2. Sensing volume of the cosmic-ray probe comprises a hemisphere in air (of radius R) and a 

22 cylinder in soil (of height D). All hydrogen within the sensing volume is reflected in the 

23 measured neutron intensity. The horizontal footprint, R, depends on air properties: mainly 

24 density and water vapor content. The vertical footprint depends on soil properties: mainly bulk 

25 density and total hydrogen content (pore water, lattice water, organic matter water). 
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26 

27 



29 Fig. 3. Geometry of a multipath signal, for antenna height (Ho) and satellite elevation angle (E). 

30 Black lines represent the direct signal transmitted from the satellite. The gray line is the reflected 

31 signal from the ground. The solid line represents the gain pattern of the antenna. Dashed circles 

32 indicate relative power levels of the gain pattern. (Reproduced from Larson et al., 2008) 
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34 

35 

36 

37 

38 

39 

40 


Soil Moisture Content 


u 


100 


130 


160 


WCR 

GPS 

Precipitation 



190 220 

day of year 


250 


260 


310 


Fig. 4. Soil volumetric water content (YWC) measured by five water content reflectometers at 
2.5 cm depth (grey lines), soil water content estimated by GPS-Interferometric Reflectometry 
(circles), and daily precipitation totals (bars) from a site near Marshall, CO, United States. 
(Adapted from Larson et al., 2010). 
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42 Fig. 5. Location of study site used by Striegl and Loheide (2012) (a), aerial photo of active DTS 

43 transect with three independent soil moisture monitoring stations (b), and schematic diagram of 

44 active DTS system components (c). Reproduced from Striegl and Loheide (2012). 
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46 
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47 Fig. 6. Four hour rainfall totals and DTS measured average temperature rise eight minutes after 

48 heating began for each 2-m interval along the 130-m cable transect (a), estimated soil moisture 

49 values based on the active DTS data (b), and a plot of active DTS soil moisture estimates and 

50 independent soil moistuture estimates versus cable position on 25 Oct. 2010 at 16:00 (c). 
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52 

53 


54 



57 Cesbio- Mira). 
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59 
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60 

61 Fig. 8. Monthly soil moisture product (September 2010) expressed in m m' . Note the wet 

62 patches in Argentina or the receding Intertropical Convergence Zone influence in Sahel. Where 

63 topography is too steep, RFI too important, vegetation too dense (tropical rain forest) or soils are 

64 frozen /covered by snow, the retrievals are either not attempted or not represented. 

65 

66 


5585 Guilford Rd„ Madison Wl 53711 


Page 97 of 110 


Soil Science Society of America Journal 



5585 Guilford Rd„ Madison Wl 53711 



Soil Science Society of America Journal 


Page 98 of 110 


70 



71 Biome 8 

72 Fig. 10. Nine AirMOSS flux sites covering major distribution of vegetation types in North 

73 American biomes. 

74 

75 
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77 Fig. 11. AirMOSS flight path made up of four flight lines, Metolius flux site, Cascade 

78 Mountains, Oregon. 
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82 Fig. 12. AirMOSS three band (Red = HH, Green = HV, Blue = YV where H is horizontal 

83 polarization and Y is vertical polarization) raw data image showing the spatial variation of soil 

84 moisture over the Metolius flux site, Cascade Mountains, Oregon along with soil roughness and 

85 vegetation effects which have not yet been removed. Volcanic feature in center of image is 

86 Black Butte cinder cone. 
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100 

101 Fig. 13. In situ soil moisture monitoring sites across the Continental U.S. 
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102 



103 

104 

105 Fig. 14. Overview of soil moisture stations currently contained in the International Soil Moisture 

106 Network (ISMN). Green dots show the stations that are still measuring soil moisture, red dots the 

107 stations that were imported from the Global Soil Moisture Data Bank. 
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KdppervGeiger Climate Class 

Fig. 15. Number of stations found within and area covered by the different Koppen Geiger 
classes after Peel et al. (2007). For the class legend we refer to the original publication. 
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113 

114 

115 

116 
117 



Fig. 16. Water deficit estimation by the atmospheric water deficit (AWD) method and soil water 
deficit methods for the 0- to 40- (SWD40) and 0- to 80-cm depths (SWDgo), with corresponding 
water deficit thresholds. Averages of 15 yr for Hollis, OK. (Reproduced from Torres et al., 
2013). 
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Fig. 17. Drought probabilities estimated by the AWD method and SWD methods for the 0- to 40- 
(SWD 40 ) and 0- to 80-cm depths (SWDgo). Average for 15 yr and eight sites in Oklahoma for 
days of the year 121 to 304. (Reproduced from Torres et al., 2013). 
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125 Fig. 18. Total rainfall (in), total reference short-crop ET (in), and average plant available water 


126 (in) across Oklahoma in May 2012 as measured by the Oklahoma Mesonet. 
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128 

Downward 



Negative Feedback 

◄ 

129 

130 Fig. 19. Schematic of principle atmospheric boundary layer interactions with the land surface 

131 conditions (modified from Ek and Mahrt, 1994 and Ek, 2005). Note that two consecutive 

132 negative feedbacks result in a positive feedback. 
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134 


135 

136 

137 

138 

139 

140 

141 



Fig. 20. Daily mean for each day in July 2006, averaged over Europe, of the observation minus 
6-hour forecast of a) screen-level temperature (K), b) screen-level relative humidity (%), and c) 
near-surface soil moisture (m 3 m‘ 3 ), from i) no assimilation (black, solid), and assimilation of ii) 
screen-level temperature and relative humidity (black, dashed), iii) AMSR-E near-surface soil 
moisture (grey, solid), and iv) both (grey, dashed) experiments. The assimilation was performed 
with an EKF using Meteo -France's ISBA land surface model. 
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144 Fig. 21. Contrasting soil water depletion profiles from Central Site East-Tiger Bush, HAPEX- 

145 Sahel project a) vegetated section and b) bare soil section (modified from Cuenca et al., 1996). 


146 
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149 Fig. 22. Cumulative runoff for Niccone at Migianella catchment: observed (Qobs), simulated 

150 (Qsim), and simulated with assimilation of ASCAT soil wetness index (Qass). (Reproduced from 

151 Brocca et al., 2010). 
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